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ABSTRACT 

Heat exchangers with one shell pass and n tube passes 
are often referred to as 1-n exchangers. The heat transfer 
literature contains many references to studies of 1-n 
exchangers when n is even but apparently little work has 
been done with respect to the 1-n exchanger when n is odd. 
This thesis greatly expands the theoretical study of 1-n 
exchangers with n being odd. While a completely closed form 
solution was found to be unfeasible, a polynomial approxi- 
mation has been developed that yields the effectiveness ( e) 
of the two possible arrangements of the 1-3 exchanger as a 
function of the capacity rate ratio (R) and the number of 
transfer units (N-^u). It is also shown that the effective- 
ness of the arrangement with two counterflow and one 
parallel flow tube side passes exceeds that of some of the 
1-n exchangers with n even. 
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NOMENCLATURE 



English Letter Symbols 

A = Exchanger heat-transfer surface, sq m 

A m = Coefficient of the mth value dimensionless 

Ao = Coefficient, dimensionless 

Ai = 1 st order coefficient to be multiplied by N-^ u , 
dimensionless 

A 2 = 2 nd order coefficient to be multiplied by N tu 2 , 
dimensionless 

A 3 = 3rd order coefficient to be multiplied by N^u^, 
dimensionless 

A 4 = 4th order coefficient to be multiplied by N tu 4 , 
dimensionless 

A 5 = 5th order coefficient to be multiplied by N tu 5 t 
dimensionless 

a = Exchanger heat-transfer surface, sq m/m 

C = Capacity rate, W°/K. Also designates dimensionless 
arbitrary constant 

Cp c = Specific heat at constant pressure of cold fluid, 
J/kg°K 

C D h = Specific heat at constant pressure of hot fluid, 

J /kg° K 

D = Empirical value of effectiveness (computer generated), 
dimensionless 

e = Error. Also used as the exponential function 

F = Logarithmic mean temperature difference correction 
factor, dimensionless 

L = Exchanger length, m 

N = Number of effectiveness empirical data points used to 
determine a curve for R, dimensionless 
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n = 

n c = 
n h = 
N t u = 
P 

q 

Qmax = 

R 

S 

Sr = 
T 

T pi = 

Ti = 

t 2 = 
tl = 
t2 = 
*a = 
t ab = 
tb = 
t bc = 
t c = 
U 

w 

x = 



Number of tube passes, dimensionless. Also number of 
equations, dimensionless 

A related Nt u per unit length hot side, m”l 
A related Nt u per unit length cold side, m"l 
Number of transfer units, dimensionless 
Temperature group, dimensionless 
Total rate of heat transfer, W 
Maximum total rate of heat transfer, W 
Capacity rate ratio, dimensionless 
Temperature group, dimensionless 

Sum of the squares of the residuals, dimensionless 

Hot fluid temperature, °C 

Particular integral, dimensionless 

Hot fluid temperature in, °C 

Hot fluid temperature out, °C 

Cold fluid temperature in, °C 

Cold fluid temperature out, °C 

Cold fluid temperature 1st pass, °C 

Cold fluid temperature between 1st and 2nd passes 

Cold fluid temperature 2nd pass, °C 

Cold fluid temperature between 2nd and 3rd passes 

Cold fluid temperature 3rd pass, °C 

Overall heat transfer coefficient, W/m^ - °C 

Mass flow, kg/sec. Also the product of w and L, 
dimensionless 

length coordinate, m. Also used to represent 
a constant value in a sequence 
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y = Sum of mth degree polynomial, defined by eq . (53), 
dimensionless 

Z = A product of z and L, dimensionless 
z = A related N-^u P er unit length, hot side, 1/m 

Greek Letter Symbols 

a = Root of auxiliary differential equation, 1/m 

e = Exchanger effectiveness, dimensionless 

X = Combination of variables defined by equation (11), 
dimensionless 

a) = A related N-tu per unit length, hot side, 1/m 

<j> = A combination of terms defined by eq. (38), 

dimensionless 

E = Summation, dimensionless 

a 2 = Variance, dimensionless 

6 m = Mean temperature difference for exchanger, °C 
9 = Indicates partial derivative, dimensionless 

Subscripts 

c = Cold fluid 

h = Hot fluid 

i,j,k = Values in a sequence 

m = Degree or order, an exponent 

1 = inlet 

2 = outlet 

Special Symbols 

[A] = An m x n matrix, symmetric 
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[K] = An m x 
f L ] = A lower 
[T] = An m x 
[ 0 ] = A nu 1 1 



n matrix, symmetric 
triangular matrix 
1 vector 
vector 
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INTRODUCTION 



I . 

A. BACKGROUND 

When analyzing the standard counterflow heat exchanger, 
it becomes apparent that, from a practical standpoint, it is 
often difficult to obtain a high velocity for one of the 
fluids when this fluid is constrained to flow through all of 
the tubes in a single pass. This leads to a possibility of 
a low overall heat transfer coefficient which cancels the 
advantage of the high logarithmic mean temperature 
difference which is obtainable in true counterflow. 

The quest for flow arrangements for increased heat 
recovery has led to arrangements that yield increased tube- 
side velocities and higher overall heat transfer 
coefficients even at the expense of a departure from the 
ideal true counterflow arrangement. Thus, the design may be 
modified so that the tube side fluid is carried through 
fractions of the tubes consecutively. 

Heat exchangers of this type with one shell pass and n 
tube passes are often referred to as 1-n exchangers. These 
exchangers, such as the one shell pass, two-tube pass (1-2) 
parallel-counterflow exchanger (see Figure 1.1), are 
configured such that all of the tube side fluid flows 
through the two halves of the the tubes successively. A 
single channel is employed with a partition to permit the 
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Figure 1.1 1-2 Parallel-Counterflow Exchanger 



entry and exit of the tube side fluid from the same channel. 
Note the baffles used to induce turbulence causing the 
liquid to flow through the shell at right angles to the axes 
of the tubes thus helping to create a higher shell side 
velocity with higher shell side heat-transfer coefficients. 

To date, much work has been done on finding the true 
logarithmic mean temperature difference for heat exchangers 
with an even number of tube passes. Little work, however, 
has been done with regard to heat exchangers having an odd 
number of passes. This is primarily due to the method 
employed in deriving an analytical solution to measure the 
overall effectiveness of a heat exchanger. 

Exhangers with an even number of tube passes often 
present a configuration problem especially in a marine appli- 
cation where the inlet and outlet of the cooling fluid must 
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be one the same side of the exchanger header (see Figure 
1.1). This particular problem could be alleviated by going 
to heat exchangers with an odd number of tube passes. 
Currently, precise mathematical expressions for the 
effectiveness of the 1-3 exchangers do not exist. Hence, a 
theoretical examination is reported on here which considers 
the effectiveness of the 1-3 parallel-counter flow exchanger 
which is shown in Figure 1.2. 

Before doing this it is important to mention the work 
that lead to the effectiveness method and the development of 
work on heat exchangers with an odd number of tube passes. 




Figure 1.2 1-3 One Shell Pass Three Tube Pass Exchanger 
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This will be considered in Section II. Kern [Ref. 1: 
pp. 224-226] , makes the interesting point that the optimum 
exchanger requires an exchanger capable of providing the 
optimum fluid-flow velocities on the shell as well as the 
tube sides. This might frequently entail the use of an odd 
number of tube passes or an odd tube length. 



B. WHY EFFECTIVENESS AS A FUNCTION OF N tu 

It is also noted that Kays and London [Ref. 2: pp. 

24-29] indicated that the effectiveness as a function of 
^tu ( e ~ Nf U ) method is the favored approach for evaluating 
a heat exchanger's performance because: 

1. The effectiveness value stands alone as a dependent 
variable and should not appear directly in the 
abscissa and indirectly in the ordinate of a graphical 
display . 

2. The log-mean difference equation misleadingly 
simplifies the notion of what is involved in heat 
exchanger design theory, since the implication is that 
only a rate equation is required. 

3. The e - N tu approach simplifies the algebra involved 
in predicting the performance of complex flow 
arrangements . 

4. The more meaningful arguments are related to ease of 
use in design work. Two prime examples of these are: 

a) Given the overal heat transfer coefficient, U, the 
two fluid capacity rates, C c and C^ , and the 
terminal temperatures, determine the required 
surface area, A. 

b) Given A, U, C c , C^ and the inlet temperatures of 
both streams, determine the outlet temperatures. 
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THE DEVELOPMENT OF THE EFFECTIVENESS METHOD 



II . 

A. LITERATURE SURVEY 

It was Nagle [Ref. 3: pp . 604-609], in 1931, who 

credited Davis [Ref. 4] with a simplified method for 
computing actual temperature differences between two heat- 
interchanging streams which depart from true counter or 
concurrent (parallel) flow. This is now the familiar "F 
factor" method which expresses the actual mean temperature 
difference 0 m in q = UA 9 m as a fraction F of the counterflow 
logarithmic mean temperature difference, LTMD, 0,^ via 9^ = 

F ®mc • 

The example of initial interest was the 1-2 exchanger 
with a single shell pass and two continuous tube passes in 
counter and concurrent flow with it. The method involved 
derivation of the actual temperature difference for the flow 
pattern and formed the ratio F = 9m/ ^mc • This familiar LMTD 
correction factor was plotted conveniently as functions of 
the effectiveness, e, and the capacity rate ratio P with R 
as a parameter. These mean temperature difference 
correction charts are available for many flow arrangements 
[Ref. 1: pp . 829-833 and Ref. 5]. The effectiveness, e 

(often called P or S) , is always the cold fluid 
effectiveness and R is always the capacity rate ratio of 
cold fluid to hot fluid. 
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Nagle detailed assumptions and derivations for the 1-2, 
1-4 and 1-6 exchangers. The F factors were obtained by 
Nagle through graphical integration and were accompanied by 
the comment that F factors for the 1-2 exchanger could be 
applied with negligible error to 1-4 and 1-6 exchangers. 
Underwood [Ref. 6: pp. 145-148] rederived the equations of 
Nagle for 1-2 and 1-4 exchangers to eliminate the need for 
obtaining F factors by graphical integration. 

Bowman [Ref. 7: pp. 541-544 ] pointed out that for a 
very large or infinite number of tube passes, the F factor 
approached, as a limit, its value in crossflow with both 
fluids completely mixed. It was further stated that even at 
the limit, the F factors were only 1 to 2 percent lower than 
those for the 1-2 exchanger. A previous paper by Kraus and 
Kern [Ref. 8] did not confirm the generalization that 1-n 
exchangers differed only negligibly from the 1-2 exchanger 
although this lack of confirmation was obtained on an 
e = f(R, Nt u ) basis. Moreover, the Kraus-Kern work does not 
confirm the generalizations on an F = f(R, N^u , e) basis. 

From the standpoint of usefulness and good accuracy, it 
is essential that F factors, if they are to be used in 
preference to e = f(R, N^u» flow arrangement), be obtained 
with precision. Plots of F = f(R, £, = P or S) [Ref. 1: 
pp. 829-833 and Ref. 5] show that the curves for particular 
values of R approach infinite slope as F decreases. While 
this can be partially alleviated by restricting R. < 1.0 (a 
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constraint used in the e = f(R, N-(- u , flow arrangement) 
approach), it is seen that small errors in the interpolation 
for R or e = P or S can result in large fluctuations in the 
value of F. 

In a comprehensive paper, Bowman, Mueller and Nagle 
[Ref. 9: pp. 283-294] presented graphs of F factors for 
shells with one through six shell passes and numbers of 
continuous tube passes respectively double the number of 
shell passes. In view of the earlier references to Nagle 
and Bowman, it should be noted that F factors were computed 
for the 1-2 exchanger in [Ref. 9: pp . 283-294] using the 
equations of Underwood [Ref. 6: pp. 145-158]. 

Ten Broeck [Ref. 10: pp. 1041-10421 prepared a graph 
of the dimensionless groups now known as e, R and N^ u for 
the 1-2 exchanger. Such a graph had the added versatility 
of simplifying the calculation of performance in a given 
exchanger when operating at conditions different for those 
for which it was designed. Kays and London [Ref. 2: 
pp. 63-74] prepared similar graphs and tables of 
e = f(R, Ntu> flow arrangement) for the 1-2 exchanger and 
for several cases of crossflow and periodic flow. 

The foregoing describes the early history of the search 
for the so-called Logarithmic Mean Temperature Difference 
Correction Factor, F, with regard to heat exchangers having 
an even number of tube passes. It is a fact, however, that 
certain space economies could be realized from exchangers 
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having an odd number of tube passes so that the tube side 
fluid could enter and leave the exchanger at opposite ends 
of the exchanger (see Figure 1.2). 

B. FISCHER'S WORK 

With the foregoing in mind, an extensive search has been 
conducted to obtain e - N-^u data for the so called "1-3" and 
"1-5" exchangers. This search has uncovered a single work, 
that of Fischer fRef. 11: pp. 377-383], which summarizes 
the historical development covered here and contains only a 
small section on the 1-3 exchanger. This work by Fischer 
develops an equation for true mean temperature difference of 
the 1-3 exchanger and casts the results in terms of F rather 
than e. Moreover, the work treats only the case where the 
three tube passes are arranged with two in counterflow and 
one in parallel flow (1-3:2C) making no mention of the one 
counterflow and two parallel flow (1-3:2P) case (see Figures 
2.1 and 2.2). In addition, the equation developed to yield 
F must be solved using a trial and error solution. 

The present work is aimed at continuing the Fischer 
investigation for several reasons: 

1. A solution is needed for effectiveness, e, as a 
function of capacity rate (R) and number of transfer 
units (N tu ). 

2. This solution should be in a closed form if at 
all possible so that it will be computationally 
efficient and useful in both the design and analysis 
frameworks . 
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3. Because valid design data can evolve from a 
polynomial approximation. The search should not be 
abandoned just because a closed form solution does not 
result from an analytical approach. 

4. Data is needed for the (1-3:2P) two parallel-one 
counterflow configuration. 

5. A 1-3 exchanger in a marine (shipboard) application 
may result in a considerable space saving over its 
1-n counterpart with n even. This would be evident 
on the outside of the exchanger where it would be 
immediately noted that the 1-3 exchanger has tube 
side inlet and outlet at opposite ends of the 
exchanger . 

The next section confirms Fischer's result and shows 
that a closed form solution cannot be obtained for the 
effectiveness of the 1-3 exchanger. Sections IV and V 
demonstate how, through numerical analysis assisted by a 
computer, a polynomial solution can be derived that will 
yield the effectiveness to engineering accuracy. 
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Ill . 



AN ATTEMPT AT A CLOSED FORM SOLUTION 



A. EFFECTIVENESS AS A FUNCTION OF CAPACITY RATES AND 

EXCHANGER SIZE 

This section deals with an investigation into the 
effectiveness, e, of a one shell pass and three tube pass 
heat exchanger, whereby e compares the actual heat transfer 
rate to the thermodynamically limited, maximum possible heat 
transfer rate as would be realized only in a counter flow 
heat exchanger of infinite transfer area. This exchanger 
heat transfer effectiveness is given by 



e 



_q 

q 

^max 



C h^ T hot,in T hot,out^ 

C TT ! - t ) 

min v hot, in cold, in' 



^c^cold ,out t cold,in^ 

C CT - t 

min v hot, in cold, in) 



where C m i n is the smaller of the C^ and C c magnitudes. 

Thus, e possesses the significance of effectiveness of the 
heat exchanger from a thermodynamic point of view, with the 
magnitude of the effectiveness completely defining the heat 
transfer performance. In general we express e = f(N-t u , R, 
and flow arrangement) and when the flow arrangement is under- 
stood, it is said that e = f(N-t u , R) . [Ref. 2: pp . 14-26]. 

The number of heat transfer units N tu is a nondimensional 
expression of the "heat transfer size" of the exchanger. 

When N-£ U is small the exchanger effectiveness is low, and 
when Nt u is large, e approaches the limit imposed by the flow 
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arrangement and thermodynamic conditions asymptotically. 



From inspection of the definition of N-t u 



A 




C . 
min 



1 



S DdA 

o 



it is clear that the overall conductance and transfer 
area affect the costs of attaining a high value for N-t u , 
ergo high e. The capacity rate ratio, R, as defined by 



is simply the ratio of mass flow rate times specific heat 
capacity for the two streams. These can be considered as 
flow stream thermal-capacity rates, i.e., energy storage 
rate in the stream per unit of temperature change. [Ref. 2: 
pp. 14-26] 

The attempt taken in this thesis to develop a closed 
form solution has used the basic fundamentals of heat 
transfer as well as those indicated above. A closed form 
solution for e was sought for both 1-3 exchangers with one 
having two out of three tube passes in parallel flow and the 
other having two out of three tube passes in counterflow. 

The analytical approach taken, and demonstrated in this 
section, is for two out of three tube passes in counterflow. 




C 

min 



max 
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B. ANALYTICAL DEVELOPMENT 



The derivation for the effectiveness, e, of the 1-3 
exchanger as a function of the capacity rate ratio, R, and 
number of transfer units, N-t u , depends on several 
assumptions . 

(1) The overall coefficient of heat transfer, U, does 
not vary within the exchanger. 

(2) The specific heat of both hot side and cold side 
fluids does not vary. 

(3) Each fluid is thoroughly mixed, that is, the 
temperature of both hot and cold side fluids is 
uniform over any cross section. 

(4) Steady flow conditions are maintained. 

(5) Heat losses to or from the environment are 
negligible . 

(6) No change of phase takes place; all heat 
transferred is sensible heat. 

(7) There is equal heat transfer surface in each pass. 

The configuration is shown in Figure 3.1 where the three 

tube passes are designated with subscripts a, b and c. The 
temperature of the hot (shell side) fluid is indicated by 
upper case letters. For the cold (tube side) fluid, lower 
case letters are used. The subscript 1 always refers to the 
fluid inlet and the subscript 2 always refers to the fluid 
out let . 

With Wft and Cp^ designating mass flow (kg/sec) and 
specific heat ( Jou les/kg - 0 K) of hot fluid entering at T-^ and 
leaving at T 2 we define a capacity rate for the hot side 

^h = ^hCph 
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T 



1 



t 



2 



t 



ab 



x :0 



x:L 



Figure 3.1 Two Counter One Parallel Configuration for 
Development of Effectiveness Relationship 
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In similar fashion for the cold side (with W c and C pc ) 
entering at ti and leaving at t2, we have 

C c ~ W c C pc 

We then obtain an energy balance for the entire 
exchanger 

c h( T l ~ t 2) = C c (t2 ~ ti) (1) 

Over the right hand side of the exchanger (Figure 3.1) 
^h( T l — T) = C c (t2 ” t c + t^-, — t a ) (2) 

and a differentiation gives 

C h dT = C c (dt c - dtfe + dt a ) (3) 

Across dx , with a (m^/m) , the surface per running meter 
of length of pass so that A = 3aL is the total surface in 
the exchanger, we may write the heat transferred to the 
element dx in each cold pass. 





Uadx(T - 


^a) 


(4a) 


C c dt b = 


-Uadx(T - 


tb) 


(4b) 


C c dt c “ 


Uadx (T - 


tc ) 


(4c) 



Here it should be observed that due cognizance has been 
taken of the direction of the flow in each cold fluid pass 



30 



with respect to the positive sense of the length coordinate 
x, and U is the overall heat transfer coefficient (W/m2-°C) 
With eqs. (4) in eq. (3) 

C^dT = Ua(3T - t a - t^ - t c )dx 



or 



dT 



-j — = n, ( 3T - t - t, - t ) 
ax h v a b c 



(5) 



where 



n - Ua 

n h XT 
h 



is a sort of N tu per unit length for the hot side. 
Now differentiate eq. (5) 



j 2m jm dt dt K dt 

d T _ f Q dT _ a _ b _ c N 

, 2 ~ n h^ ox dx dx dx ' 

dx 



and with eqs. (4) substituted 



= 3n. ^ - n n. (T - t + t. - t ) 
,2 hdx ch v a b c 

dx 



( 6 ) 



where 



n 



c 



. Ua 

XT 

c 



where again the resemblence of n c to N^ u can be noted. 
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From eq . (2) we obtain 






(T, ~ T) ~ t„ = t. - t - t 



(7) 



and with eq . (7) put into eq . (6) we obtain 



d 2 T . dT 

~~ ! 1 ~ Jn h dx 
dx 



- "c"h > T + <T < T 1 ' T > 
c 



V 



or 



where 



- 3n. 4^ = -n n. ^ f(R - 1)T + T, - R t 0 ] (8) 

,2 hdx chC' v c ' 1 c 2 J v ' 

dx c 



R C “ ^c/^-h 

is the capacity rate ratio for the cold side 
Notice that 



C h Ua Ua C h ,Ua 2 

n c n h Tr"Tr*ir , u-- ( n- ) " m 
c c h c c 



and 



R h R 



XT 

c 



a capacity rate ratio for the hot side. Then, algebraic 
adjustment provides 
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(9) 



d 2 T 



dT 



1 - R. 



dx 



^ n h Hx* + m ^ d “ Ti ) 



h 



R h 1 



which is a linear, non-homogeneous , second order 
differential equation with constant coefficients having a 
complementary function 



a-, x oirtX 

T c = C l e + C 2 e 



( 10 ) 



where Ci and C 2 are arbitrary constants and where 



3n, 



h .1 , 0 2 . . 11 . , 

a l ’ a 2 5“ " ~2 9n h 4m( R ^ 



1 - R h 1/2 



h 



But 



3n h . n h ,n _ 4m 

2 " 2 1 2 ' 

n. 



1 - R. 1/2 



R, 



■)] 



,2 (C. C. 2 0 . 

m _ (Ua) v h' _/h. _ R 2 _ 1 

2 .2 /TT ,2 ( 'C ; h _ 2 

n. (C ) (Ua) c R_ 



so that 



a l » a 2 



3n u n. 

£ + Jl rg 

2 ~ 2 



? 1 - R h 1/2 

V 



or 



a l * a 2 



2“ 



(3 ± X) 



( 11 ) 
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where 



A = [9 - 41^(1 - VI 172 (12) 

Designate the particular integral as T pi and by the 
method of undetermined coefficients let T p i = P so that in 
eq. (9) 



1 - R. t 9 

m( — - ) P = m( p— - T ) 

h h 



This makes 



T . 
Pi 



1 0 R h 

p - 1 TC * T l' 

h h 



so that 



T . 
Pi 



t 2 R h T l 
1 - R. 



(13) 



The general solution to eq . (9) is the sum of eqs. (10) and 
(13) 



T(x) 



cx*< x oCrt x 

C l e + C 2 e 



+ 




R h T 

R 



(14) 



where the arbitary constants, Ci and C 2 are evaluated from 
conditions at x = 0 and x = L. At x = 0 , T(x =0) = T 2 and 
at x = L, T(x = L) = Ti . When these are inserted, in turn, 
into eq . (14), one obtains a pair of linear algebraic 

equations in the unknowns Ci and C 2 
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to show that 



T 2 - c i + c 2 + T pi . 



cu L ciqL 

T 1 - °l e + C 2 e + T pl 



where T p ^ is given by eq. (13). 

It is only a matter of algebra 



C, = 



(T, - T . ) - (T„ - T . )e 

v 1 pi 7 _2 pi 7 

a^L a 2^ 
e - e 



a 2 ] 



(15a) 



and 



C„ = 



(T 0 " 



T . )e 
Pi 



a^L 



- (T, - 



T . . 
Pi) 



a^L otgL 

e - e 



It is easy to see from eq . (1) that 




(t 2 t 1 ) 

( T 1 - V 



so that 

t 2 = ti + Rh(Ti - T 2 ) 

Use of this in eq . (13) shows that 




t l + R h^ T l ” T 2 ) R h T l 



1 - R, 



(15b) 
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or 



t l R h T 2 



Pi 









(16) 



indicating two alternative forms for T p i given by eqs. (13) 
and (16). 

Insertion of eqs. (13) and (16) in eqs. (15) for and 
C 2 will yield after some algebra 



, T 1 t 2, _ , T 2 t l N ^ a 2 I 

( 1 - R, ^ ( 1 - R, ) e 



C. = 



V 



h 



a^L 

e - e 



(17a) 



and 



/ T 2 t l N a l L , T 1 t 2 N 

( 1 - - R -) e " ( r — D ) 



C 2 



h 



TT - R 



h 



a^L a^TT 

e - e 



(17b) 



Equation (14) is an expression for the hot side 
temperature at any location in the exchanger in terms of the 
extreme temperatures, t\, t2, T]_ and T 2 * 

Next take eq . (5) and set it equal to the derivative 
of eq. (14) noting that Ci , C 2 and T pi are all known 
constants . 



33E = V 3T ' - ‘b - V - “l 0 !® 01 * + “2 C ! eV 08) 

At x = 0, where T = T 2 , t a = t^ and t^ = t c = t^ c 
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( 19 ) 



dT 



3 x " n h (3T 2 t l 2t bc ) " °1 C 1 + a 2 C 2 



and if we subtract eq . ( 4 a) from eq. ( 4 c) we obtain 



dt - dt 
a c 

t - t 
a c 



c 



-n dx 
c 



which can be integrated using C3 as the constant of 
integration . 



t a t c C 3 e 



-n c x 



and at x = 0 where t a = ti and t c = t^c 



tl - t bc - C 3 



or 



tbc = ti - C 3 



In addition at x = L, t a = t a b and t c = t2 so that 



t ab t 2 C 3 e 



n c L 



or 



n - t ab t 2 
u 3 r nTL 



= ^ah " t 2 )e 



n c L 



This gives a relationship between t a b and tb c 



N. 



t bc t l (t ab t 2 )e 



( 20 ) 
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where N c = n c L can be considered as the total number of 
transfer units for the cold side. 

Return now to eq. (18) and look at the conditions at 
x = L where t a = tfo = t a b , t c = t 2 and T = . These 

conditions in eq . (18) give 



where again we remember that C\ and C 2 are known constants. 
Solving for t a b 




2t 



ab 




> + 3Tj - t 2 



and with this in eq . (20) 



2t 



be 




OoL 

+ a 2 C 2 e ) + 3t 2 - 3T 1 ] + 2t 1 



Then with eq . (21) in eq. (19) 



N c 

a l C l + a 2 C 2 = 3n h f(T 2 " t i> + e ( T i “ t 2 )] 



N cl.L apL 

-e C (a 1 C 1 e + * ) (22) 



Equation (22) confirms Fischer's result fRef. 12: 



pp. 377-383] and at this point in his development he 
branches off to seek an expression for the Logarithmic Mean 



Temperature Difference Correction Factor, F. The attention 
here is focused on e = f(R, N^ u ) and the balance of this 
section continues in this vein. 

Look at the aiCi + 012 C 2 term in eq . (22). Use of eq . 

(11) allows the representation of aiCi + 012 C 2 

n. n. 

7 y— (3 + X)C 1 + ^ (3 - X)C 2 (23a) 

or 




+ C 2> 




C 2> 



(23b) 



Then from eqs. (17) 



C 1 + C 2 - 



( T 2 ~ t 1 )[e 



Ot^L 



- e ] 



a. L a^L 
(1 - R h )\e - e * ] 



or 



C 1 + C 2 - 



T 2 ~ t l 
1 - R. 



(24) 



Moreover 



C 1 - C 2 " 



ou L a^L 

2(T X - t 2 ) - (T 2 - t 1 )[e + e 1 

ol, L cuL 
(1 - R h )[e 1 - e 2 ] 



(25) 



Now let 
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(0 



(26a) 



3n 



h 



and 



Xn 



z = 



h 



(26b) 



so that 



or 



« 1 L a 2 L ( uj+z ) L (w-z)L 

e - e - e K - e v y 



goL zL ojL -zL 
= e e - e e 



wL, zL -zL N 
= e (e - e ) 



ON L Oto L t 

1 2 n ojL . , T 

e - e = 2e smh zL 



(27) 



Moreover, it is easy to see that 



OCh L 0tr>L t 

12 0 o)L , T 

+ e = 2e cosh zL 



(28) 



If eqs. (24) through (28) are collected and put into the 
expression of eq . (23b), the result is 



3n, Xn. 

a,C« + a 0 C 0 = (C 1 + 0 o ) + _J1 (C, - C„) 



‘1^1 2 2 



(23b) 



T 2 - t 1 2(T X - t 2 ) - (T 2 - t 1 )2e aL csch zL 

= « \ j = ~Z 1 + Z f 7T 1 



"h 



(1 - R^)2e U) ^ sinh zL 
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or 



a l C l + a 2°2 



^2 ~ ^ 1 
"< i . n. > + 



T - t t - t 

z [ ( — j — — — ) e ^csch zL - ( j ■ _ ) coth zL] (29a) 

n h 



and this could also be written as 



T 2 ^ 1 

a l^l + a 2^2 = ^ — j[ _ p ) [ “ “ 2 coth zL] 



V 



+ z( 



T 1 t 2 x - ojL 



1 - R, 



)e csch zL 



(29b) 



The next step is to reduce the right hand side of eq . 
(22). Use of eq . (11) permits the representation 



oil L doL 3n. in. . 3n^ Xn^ ou L 

o 1, ^ z /h, n x ~ a.L , , n h . T. 

a l^l^ a 2^2 e — — 5 — — 2 — )C^e 1 ■*" ( — g — — — 2 — 



or with eqs. (26) inserted 



a. L a^L a 1 L a~L cl. L a^L 

a l C l e + a 2 C 2 e = °^ C i e + C 2 e ) + z ( c 1 e “ c 2 e ) 



(30) 



But by eqs. (11) and (17) 



a.L oujL ( a-. + ap ) L (a, + a^)L 

a^L a 2 L (T^t^e - e Z ) + (T 2 “t 1 )[e - e ] 

C l e + C 2 e = aJT ^ L 



(1 - R h )(e - e “ ) 
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or 



CUL OtcrjL T 1 t n 

C l e + C 2 e - I^-fT 2 



(31) 



a^L a^L 

- C^e may also be simplified. Again using 



eqs . (11), (17a) and (17b) 



ou, L cv>L (cu + ar,)L ( cu + Op )L 

a x L a 2 L_ (T 1 -tg)(e + e Z )-(T 2 ~t 1 )fe + e * ] 

C l e “ C 2 e = aTE ouE 

(1 - R. )(e - e ) 



The exponential term at the far right in the numerator 
is really quite simple. From eq . (11) 



3n, An. 3n^ An. 

, n n v , n n x 0 

a l a 2 ” ( ~T~ ~2~ ) ( ~2~ ~2~ ) " 3n h 



and by eq. (26a) a\ + a 2 = 3n^ = 2 u. Thus with the 
combination of exponentials in the numerator and the 
denominator given by eqs. (27) and (28) we find that 



a 1 L a„L 

C^e 1 - C 2 e 



2 ( T 1 - t 2 )e ajL cosh zL - (T, - t 1 )2e 2aij 
2(1 - \)e“ L sinh zL 



or 

a l L a 2 L / T 1 " t 2^ T / T 2 t l x uL u j /OON 

- C 2 e = ( y _ p — ) coth zL - ( ^ - _ R — )e csch zL (32) 
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Now with eqs. (31) and (32) put into eq . (30) we 

obtain 



„ ~ a l L ^ „ ^ a 2 L _ , T 1 t 2 x 

a l C i e + a 9 C 9 e “(l — R ) 



‘2 2 



1 - R, 



T 1 t 2. 



T - t 

+ z( ( ^ x _ p - ) coth zL - ( - - ) e csch zL] 
n K h 



or 



a-. L 



1^ a 2 L Tj. - t 2 



a l^l e + a 2^2 e = ^ 1 ~ R — ) [ w + z coth zL] 

n 



T 2 t l v wL 



- z(t g — )e csch zL 

1 " *h 



(33) 



With eqs. (29b) and (33) inserted into eq. (22) we 
obtain 



T p t < T. tn _ i 

(t= g — ) [ a) - z coth zL] + z(t K — )e csch zL 

1 " R h 1 " R h 



N 



3n h f (T 2 " V + e ^ (T 1 - t 2 ) - 



—N T i tp T p — t * t 

e c f(-» b — ) ( w + z coth zL) - z(- r - ^ — 5 — ) e csch zL] 

1 " R h 1 " R h 

(34) 
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We wish to develop an expression for the exchanger 
effectiveness so we designate the hot side effectiveness as 



T 1 - T 2 
e h “ Tpr-q 



(35) 



and we begin by simplifying eq. (34) by dividing throughout 
b y ( t 2 “ ti)/(l - Rh) to obtain 



( to - z coth Z) + ze csch z( 



^1 x _2 

T 2 " t l 



N c T i ~ t 2 

) = R + e C R( t 1 _ / ) 
2 1 



N c ” 1 2 

- e c [(a) + z coth Z) (=i 

X 2 T 1 



w 

) - ze csch Z] (36) 



where 



and 



u)L 


(37a) 


zL 


(37b) 


3n h (l - R h ) 


(37c) 



We can then let 



T 1 *2 
T 2 - *1 



and get eq . (36) to look like 
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( oj - z coth Z) + (ze w csch Z) <j> = 



N -N (W-n ) 

R + (e C F) <p - fe c (cu + coth Z ) 1 — (ze c csch Z) 



It is now a matter of algebra to solve for <j> 



R + ze 

9 “ 



(W-N ) 



csch Z - ( to - z coth Z) 



ze ^ csch Z-Re C + e c (aj+z coth Z) 



N 



N 



( 38 ) 



The next step is to represent <|> as a function of ey,. 
This is done with some algebraic gymnastics as follows: 



= 



T 1 t 2 



T 1 t 2 
T 2 " h “ T 2 ' T 1 + T 1 “ *1 



T 1 t 2 

T 1 - T ? 

< T i - t 1 )r T7 _r_ t7 1 



or 



* 



( T 1 




)(i - 



'h 



) 



Moreover 



* = 



T 1 r 2 



(T x - t x )(l - e h ) 



T 1 “ t l + r l ” t 2 

(^ - V^ 1 - V 



to ^ i ^3 ~ 

(t - t ) r i - (- t - -_- F ~) c T - f -)i 

i i ‘l X 1 A 1 *2 

(T - t )(1 ~ e ) 
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But Rh = (t 2 ~ ti)/(Ti - T 2 ) so that 






1 ■ e h R h 

_T ^T 



or 




where <(> is given by eq . (38) 

The neatness of the form of eq . (39) is deceptive 

because, unfortunately, it cannot be used to determine a 
unique value for e^. The reason for this can be found in an 
inspection of eq. (38) which provides the value of 4> which 
is used in eq . (39). 

Notice in eq . (38) that Z, W and N c are all functions of 
the product aL. On the other hand, w, z and R are functions 
of a only. Thus, it is impossible to vary a and L indepen- 
dently and still achieve a unique solution. 

For example, suppose a = 50 m^/m and L = 5 m so that 
aL = 250. A value of <|> may be obtained from eq . (38) using 

these values. However, if a = 100 and L = 2.5 so that aL is 
still equal to 250, an entirely different value of <j> is 
obtained because a = 100 rather than 50. 
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One should resist the temptation to multiply numerator 
and denominator of eq . (38) by L thereby creating a 

situation where only Z, W and N c appear along with a new 
R' = RL . Such a procedure is doomed to failure because in 
dealing with an equation derived from n equations in n+1 
unknowns, one cannot create the n+1^* 1 equation by 
multiplying one of the n equations by a constant. This 
makes the n+l th equation so obtained linearly dependent on 
one of the original n equations and the entire set becomes 
linearly dependent. 

This section represents an attempt to obtain e = f(R,N-tu) 
and the attempt has not been successful. It is now time to 
turn to the computer and this will he done in Section IV. 
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IV. NUMERICAL AND COMPUTER ANALYSIS 



With a closed form solution for e as a function of R 
and N^ u not attainable as indicated in Section III, it 
becomes apparent that an alternative method is needed to 
determine the effectiveness for the 1-3 :2C and 1-3 :2P heat 
exchangers. Kern and Kraus [Ref. 12: pp. 306-360] describe 
a computer code for a thermal analyzer. This code (program) 
makes use of node equations generated by finite differences 
and it employs a Cholesky LU decomposition scheme. 

The Cholesky decomposition, as explained by Stewart 
[Ref. 13: pp. 134-144], is best used when decomposition in 
the presence of positive definite matrices is requested. 

This is the case at hand where one tries to solve 
numerically for the temperatures that lead to the 
effectiveness of the 1-3 :2C and 1-3 :2P heat exchangers. 

A. THERMAL ANALYZER TVSSI 

The computer program employed is Program TVSSI (Appendix 
A) which is an adaptation of the thermal analyzer program 
called TVSS2 and listed by Kern and Kraus [Ref. 12: pp. 
306-360]. The adaptation consisted of changing the program 
so that it could perform the computations in the SI system 
of units and be receptive to the use of a specially created 
input file. Also it should be noted that TVSS2 was written 
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to be used in conjunction with the Honeywell H-1800 computer 
system. Therefore, it had to be modified to run on the 
Naval Postgraduate School's IBM 3033 AP system. 

The program itself is a non-linear equation solver that 
determines the temperatures at a prescribed number of node- 
points or nodes from a set of node equations in almost any 
framework (i.e., network analysis, field plotting, or fluid 
flow distribution). It has certain features that make it 
primarily an equation solver for thermal analysis. These 
features include: 

1. an ability to linearize radiation terms. 

2. an ability to allow any of the coefficients in the 

node equations to vary with temperature. 

3. an ability to provide constant heat input and heat 

input as a function of temperature at any node. 

4. an ability to consider other modes of heat transfer 

that are non-linear such as boiling and natural 
convection . 

As stated earlier, the program utilizes the Cholesky 
decomposition scheme and, because of the linearization of 
the radiation terms (a feature of the program that is used 
even though radiation does not appear in this e - N tu 
study), the program is iterative. 

Cholesky' s decomposition consists of finding a lower 
triangular matrix [L] which is capable of reducing the 
original system of equations. 

mm = roi ( 40 ) 



49 



or 



mm - fQ] = ro ] <4i> 

to the unit triangular form 

r A ] r T ] - fB] = [0] (42) 

so that the sought after elements of the column vector fT] 
can be obtained by backward substitution. 

Suppose, for example, that [K] is 3x3 and assume that 
the system [K][T] = [0] has been reduced to the form 

[ A ] [ T ] - [B] = [0]. In this event a premultiplication by 
[L] will return the system to its original form, that is 

[L]([A][T] - [B]) = [K ] [ T ] - [Q] = [0] 

This implies that 

[L][A] = [K] (43) 

and 

f L ] r R ] = ro] (44) 

These equations allow the determination of fL], [A], and 

[ B ] in a very simple manner and the matrices are uniquely 
determined because f K ] and fO] are known, or, at least are 
known after each iteration because the elements of [K] are 
linearized. For a 3x3 system 
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t A,B ] 



[K,Q] [L] 



k ll 


k 12 


k 13 


Q1 




”lll 


0 


0 




1 


a 12 


a 13 


b l 


k 21 


k 22 


k 23 


<12 


= 


l 21 


1 22 


0 




0 


1 


a 23 


b 2 


k 31 


k 32 


k 33 


<13 




1 31 


1 32 


1 33 




0 


0 


1 


b 3_ 



one may obtain the following for the elements of [L], [A], 

and [B ] . 



kil = ( 1 ) 1 i 1 + (0) l i2 + (0)l i3 = In (45) 

which shows that the first column of fL] is identical to the 
first column of [K]. 

k lj = ( 1 11 ) a l j + (°) a 2j + (0)a 3 j = luaij = knaij (46) 

which shows that the first row of fA] is equal to the first 
row of f K ] divided by r k 1 1 ] and then 

k 22 = 1 21 a 12 + 1 22( 1 ) » l 22 = k 22 “ l21 a 12 

k 23 = 1 21 a 13 + 1 22 a 23 * a 23 = ( k 23 ~ 1 21 a 13)/ 1 22 

k 32 = 1 31 a 12 + 1-32 ( 1 ) » 1 32 = k 32 “ !31 a 12 

Q2 = l21 b l + l22 b 2 » b 2 = («12 ” 1 21 h l)/ 1 22 

In the foregoing manner the elements of fL], [A] and [B] 
are obtained successively in terms of previously determined 
elements in a progression that goes horizontally from I22 
on. Thus the general relationship are seen to be 
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j-1 




Z 1 ir k r.j 

r=l 



( 47 ) 



with 



Hi = k il 



(48) 



and 




i-1 



1 . a . ] 
ir rj 1 



r=l 



(49) 



with 



a 



i 3 




(50) 



Moreover, it is observed that if [K] is symmetrical which it 
must be in our coupled set of equations (kjj = kj^) then 

1 i j ~ 3- j ilii (i»j = 1 1 2 3 , . . . n~ 1 j i ^ j ) (51) 

The modification for computation in the SI system was 
quick and simple. It involved changing a numeral in two 
places (460 to 273) and some format statements (°F to °C and 
Btu/hr to Watts). 

The conversion of TVSS2 to TVSSI for running on the 
IBM 3033 AP system (which uses the FORTVS compiler) required 



52 



a modification to the Fortran program language used in TVSS2 
in order to compile the FORTVS (basically international 
Fortran 77) system used on the IBM 3033. 

B. INITIAL MODELING 

Initially, the model was designed to find the tempera- 
tures that would allow computation of a single effectiveness 
value after a detailed set of capacity rate, coefficients 
and surface data were entered into the program. From the 
scope of the problem it was realized that multiple runs of 
the thermal analyzer program (TVSSI) would be needed. It 
therefore became necessary to develop a program that given 
Ch> C c , U, A, Ti , and ti , an input file would be created for 
use by the modified version of the thermal analyzer 
(TVSSI) . 

The first step taken was to develop a program to create 
an input file for TVSSI that would yield the effectiveness 
for a 1-4 heat exchanger which could be compared to the 
existing analytical solution for the effectiveness of a 1-4 
exchanger. With this accomplished and confidence 
established, similar programs for the 1-3:2C and 1-3:2P 
exchangers could be developed. This program was called 
NTU14 (See Appendix B) and the following parameters were 
used for all runs. 

1. 250 nodes were used. 

2. The initial temperature for the computer to begin the 
iterative process was set at 200°C. 
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3. An eventual accuracy of .05 between the final and next 
to last iterations was used. 

4. A radiation coefficient convergence factor of 0.66667 
between iterations was used. 

5. The maximum number of iterations that the computer 
was allowed to perform was set at 12. 

6. A damping factor of .8 was set as an initial damping 
based on the number of non-linear terms in all of the 
node equations. 

When the values of C^, C c , U, A, and ti are set, N^ u 
and R are then compiled and an input file for TVSSI was 
generated. In this file all node equations and internode 
conductance values were determined. This program determined 
and specified the nodes that interact with each other and 
the methods by which the interaction takes place such as 
conduction, forced convection, and fluid flow. 

The program, NTU14, makes use of the fact that each term 
in a node equation shows three things. The first is the 
node that is coupled for heat flow with the node in 
question. The second is the method of heat flow between the 
nodes. In this case forced convection and fluid flow are 
used. Finally, the node equation shows the magnitude of the 
internode heat flow. Here all the pieces of information are 
collected and presented for use by TVSSI as an input file 
with all items in the proper format. 

A comparison of the effectiveness for the 1-4 exchanger 
developed by the computer to that of using the closed form 
analytical solution for effectiveness developed by Kraus 
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and Kern [Ref. 8] as shown by equation (52) 



e = 



1 + R + 



2 [1 



4R 2 ] 1 / 2 coth( 



¥ 



tu [ 1 + 



4R 2 ] 



IT 



( 52 ) 



) + 4 tanh 



tu 



was then undertaken. The results of this comparison showed 
that over the entire range of R from .01 to 1.0 for varying 
values of N-^u from 0 to 3.25 less than a 0.5% difference was 
ever realized. A small sample of these results are provided 
in Table 1. The conclusion to be drawn here, is that the 
methodology used to develop the computer program NTU14 for 
input to TVSSI for finding effectiveness was sound and could 
then be used in the development of the 1-3 :2C and 1-3 :2P 
exchanger methodology. 



C. DEVELOPED MODELS FOR 1-3 :2C AND 1-3 :2P HEAT EXCHANGERS 
The same technique used in developing the program NTU14 
was used to generate computer programs NTU32C and NTU32P. 
These are listed in Appendices C and D. The departure for 
each of these programs from the NTU14 program is in the 
number of nodes; they are based on 200 node models as shown 
in Figures 4.1 and 4.2. An example of the output file 
generated from one of these programs is found at Appendix E 
It is these values shown in Appendix E that are used by the 
thermal analyzer to determine the temperatures T 2 and tq 
for the specific set of given initial parameters , C c , R, 
A, Ti , and ti . 
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TABLE 1 

COMPUTER TO ANALYTICAL COMPARISON 
FOR R = .5 



N tu 


ANALYTICAL 

RESULTS 


COMPUTER 

RESULTS 


% 

DIFFERENCE 


0.05 


.0482 


.0482 


0.00 


0.25 


.2094 


.2090 


0.20 


0.50 


.3569 


.3559 


0.28 


0.75 


.4628 


.4612 


0.35 


1.00 


.5398 


.5377 


0.39 


1.25 


.5963 


.5940 


0.39 


1.50 


.6379 


.6354 


0.40 


2.00 


. 6915 


.6886 


0.42 


2.50 


.7206 


.7177 


0.40 


3.00 


.7360 


.7333 


0.37 


3.25 


.7406 


.7381 


0.34 
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D. SCOPE OF COMPUTER ANALYSIS 



At this point it is possible to let the computer solve 
for temperature values that yield a value for effectiveness 
based on a specific set of initial parameters. However, 
it must be realized that many computer runs are required to 
generate enough data to ensure confidence in the results 
which cover a wide range of capacity rate ratios and N-eu 
values . 

To efficiently expedite the computer task, the Multiple 
Virtual System (MVS) with Job Entry Subsystem and Networking 
(JES3) was utilized. The MVS coupled with JES3 is more 
commonly referred to as batch processing. Based on trial 
and error, it was determined that in order to build a solid 
data base, eleven different values for effectiveness were 
needed to best represent a particular value of R. This 
was required over a range of R from R = 0.1 to 1.0 in 
increments of .01. In all, 200 curves for both the 1-3 :2C 
and 1-3 :2P exchangers were needed. This means that 2,200 
unique values of effectiveness needed to be found, plus 
100 values of effectiveness to be used for comparison with 
the 1-4 exchangers. 

To complete this task, TVSSI was slightly modified in 
accordance with the appropriate guidelines of the /job 
control language (JCL) needed to run on the batch processing 
system. These modifications are few and were needed only 
at the beginning and end of TVSSI. The modified version 
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of TVSSI has been called TVCOUNT with changes shown in 
Appendix F. It was TVCOUNT that was then used to activate 
TVSSI . 

It also became necessary to modify the three input file 
programs NTU14, NTU32P and NTU32C such that they needed to 
be compiled only once. They were then loaded in a library 
file to be used when called by another program. New 
programs utilizing the batch system were written that could 
easily be loaded with the appropriate input data for a 
specific R value. These, which are referred to as "sister 
programs," were used to go from the library file to TVSSI 
and cause TVSSI to be executed eleven times under TVCOUNT 
covering the desired range of N-t u for a specific R value. 

The revised input files called NTU14C , NTU32CC , NTU32PC and 
their associated "sister execution programs," NTU14L, 
NTU32CL, NTU32PL, are found in Appendices G through L. 

The overall system flow chart of how all of the fore- 
going is accomplished is found in Figures 4.3, 4.4 and 4.5. 
It is noted from these figures that TVSSI is referred to as 
TVSSIA through TVSSIV. These are the same programs as TVSSI 
but for bookkeeping purposes by the computer they are 
labeled A through V. 

E. COMPUTER RESULTS 

Upon completion of all data collection from the computer 
output, plots of effectiveness vs. N^ u for the whole range 
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of R were plotted and are shown in Appendices M and N. 
Thirty-three different plotting programs utilizing the 
"Display Integrated Software System and Plotting Language 
(DISSPLA)," were written to graph the data obtained. 

An example of one of these programs is provided in 
Appendix 0. 

Examination of the graphs in Figures 4.6, 4.7, and 4.8 
shows that the 1-3 :2C exchangers outperform the 1-3 :2P, 1-2 
and 1-4 exchangers. This is true in all cases, and, because 
of this, only a sample of the data was chosen to be shown in 
these figures. Furthermore, it is noted that at higher N^ u 
values, the effectiveness of the 1-3 :2C exchanger is better 
than all of the others considered, while at higher capacity 
rate ratios, the effectiveness of 1-3 :2C exchanger even 
begins to outperform the others at lower N^u values. This 
increase in performance is easily understood because it has 
been proven by Kern [Ref. 1: pp. 139-137] and others that 
greater temperature differences result when process streams 
are in counterflow than parallel flow. Thus, when there is 
a combination of the two phenomena (counterflow and parallel 
flow) occurring, then it becomes apparent that the extra 
counterflow pass must increase the exchanger's overall 
performance . 

As shown in the graphs in Appendices M and N, effective- 
ness increases as both R and N^ u increase. These curves can 
be used to give a graphical approximation of effectiveness 
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if that is all that is required. From the empirical data 
used to develop these curves, further investigation into the 
development of equations for these curves which may be used 
to determine an exact value of effectiveness when R and N-t u 
are known, is undertaken in Section V. 
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-3:2C NODE MODEL 





Figure 4 . 1 



Nodal Arrangement for the 



1-3 :2C Exchanger 
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-3:2P NODE MODEL 






Figure 4 . 2 



Nodal Arrangement for the 



1-3 :2P Exchanger 
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11 sets of initial parameters 
given (N^m. C^, C c , U, A, Tj and tj ) 
make up the input data to be 
solved for 




Figure 4.3 Computer Systems Flow Chart for 1-4 Exchanger 
Model 
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11 sets of initial parameters 
given (N tu , C h , C c , U. A, Ti and tj) 
make up the input data to be 
solved for 




Figure 4.4 Computer Systems Flow Chart for 1-3 :2C Exchanger 
Model 
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11 sets of initial parameters 
9 1 ven (N tu , C h » C c , U, A, Tj and t, 
make up the input data to be 
solved for 



rvssiL 



IE 



rvssiM 



:e 




-Generates input file 
under Batch Processing 
-Stores file in a Library 



'"-Uses Batch Processing 
-Calls input data 
-Calls Library input file 
] -Calls TVSSI 11 different 
times and designates 
as TVSSIL thru TVSSIV 



TVSSIU 



rvssiv 



T 



-Modified thermal analyzer 
-Uses Batch Processing 
-Executes TVSSIA thru TVSSIK 



-11 different corresponding 
values for effectiveness 
based on initial input data 
is listed here 



Figure 4.5 



Computer Systems Flow Chart for 
Model 



1-3 :2P Exchanger 
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Figure 4.6 Comparison of Analytical (1—2 and 1-4) to Computer 
Results ( 1-3 : 2C and 1-3:2P) at R = .2 
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Figure 4.7 Comparison of Analytical (1-2 and 1-4) to Computer 
Results (1-3 :2C and 1-3 :2P) at R = .5 
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Figure 4.8 Comparison of Analytical (1-2 and 1-4) to Computer 
Results (1-3:2C and 1-3 :2P) at R = 1.0 
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V. POLYNOMIAL REGRESSION 



A. DEVELOPMENT OF POLYNOMIAL EQUATIONS 

The empirical data obtained for the two 1-3 heat 
exchangers was designed to cover an extensive range of R 
values varying from 0.01 to 1.0 in increments of 0.01. As 
discussed earlier in Section IV, the data obtained at a 
specific value of R is the computer evaluated result of the 
effectiveness, for an associated Nt u value. With this 
accomplished, it then becomes possible to graph separate 
curves for each of the different R values as shown in 
Appendices M and N. Through a polynomial regression 
technique, as discussed in this section, it is also possible 
to develop implicit equations for the curves with e = f(Nt u , 
R). It is also apparent from an inspection of the graphical 
representation of the empirical data in Appendices M and N, 
that the curves conform to a high degree polynomial. 

However, further analytical investigation is needed to 
ascertain the exact order of the polynomial terms. This 
investigation will not only lead to the order of the 
polynomial, but to the specific equation for each curve. 

By use of polynomial regression, the least-squares 
method can be readily extended to best fit the data to the 
mth-degree for the polynomial 

y = Aq + Aj_x + A 2 X^ + ... A^™ (53) 
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with the error defined by 



ei = D i - Yi = D i - A o - Aix - A 2Xi 2 - ... - A m xi m 

where represents the empirical data value corresponding 
to xi , xi being free of error. 

The objective is to minimize the sum of the squares of 
the residuals, S r , 

m p m p 2 

S r = Z = Z (D - A - A x x x - a 2 x 2 + ... Ax; (54) 

i=l 

Because at a minimum, the partial derivatives 9S r /9Ao, 
9S r /9 A i ... 9S r / 9A m vanish, after taking the derivative of 
S r with respect to each of the coefficients of the poly- 
nomial, it can be seen that 



9S r 2 

-t-t— = 0 = -2 Z (D. - A - A,x . - A x . 
9A q v i 0 lx 2 l 



. - A x . m ) 
mi 



9S „ 

= 0 = -2 Z x . (D . - A n - A, x . - A 0 x. - ... - A x. m ) 

9A, i v i 0 li 2i mi' 



3S i 
9 A, 



= 0 = - 



2 Z x. (D. 

l v l 



A 0 A l x i 



A 2 X i‘ 



A 

-Ax. ) 
m x ' 



9S 

vr = 0 = ' 2 J x i” (D i ' A o - Vi ‘ Vi - 

m 



. - A x. m ) 
m l 
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Then by dividing by -2 and rearranging we obtain 



AqM + Ai Z + A^x ^ 2 + ... + A m X Xj[ m = Z 

A 0 Z Xi- + Ai Z Xi 2 + A2 Z Xi 3 + ... + A m X Xi m+1 = X x±r>i 

A 0 X Xi 2 + Aj_ X Xi 3 + A2 X Xi 4 + ... + A m X Xi m+2 = X Xi 2 D-L 



Aq Z Xi m + Ai X x i in+1 + A 2 X Xi m+2 + ... + A m X Xi 2m = X Xi m D i 

where all summations are from i=l through n. All of the 
foregoing m+1 equations are linear and have m+1 unknowns: 

Ao, Ai , A 2 , ... A m . The coefficients of the unknowns can be 

calculated directly from the observed data. Thus, the pro- 
blem of determining a least-squares polynomial of degree m 
is equivalent to solving a system of m+1 simultaneous linear 
equations. Putting the equations in matrix form yields 



N X Xi X Xi 2 X x ^ 3 ... X Xi m 




I 

•H 

c 


X XjL X Xi 2 X Xi 3 X Xi 4 ... X Xj L in+1 




E x ^ D 


X Xi 2 X Xi 3 X xi 4 X Xi 5 ... X Xi m+2 

• i • • • 

f • • • • 

• * • • • 


[A] = 


X XjL 2 D-L 

• 

• 


X Xi m X XjL m+ ^ X Xi m+2 X Xi m+3 ... X Xi 2m 




X Xi m D^ 



[Ref. 15: pp . 302-309 and Ref. 16: 468-474]. 
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From this point on, one finds that it is best to use a 
computer to assist in solving the simultaneous equations and 
this will also help alleviate any ill-conditioning that may 
otherwise occur. An existing "curvefit” program available 
through NON-IMSL [Ref. 16] and found in Appendix P was 
used although some modifications were made to the original 
program to best accommodate the goals of this work. 

To determine the order of polynomial that should 
eventually be used, one increases the degree of the 
approximating polynomial as long as there is a statistically 
significant decrease in the variance o^, which is computed 

by 



2 

a 



2 

Z e. 

_ i 

N - m - 1 



(55) 



In otherwords, the selection of the optimum degree poly- 
nomial is contingent upon a decreasing variance and once the 
variance begins to increase, the degree of the polynomial 
becomes too high. For all cases, it was found that the 
£ " N tu developed curves are of the 5th order. 

As shown in Figures 5.1 and 5.2 the computed values of 
effectiveness vs. N-^u for R = 0.1, 0.5 and 1.0 for both 
flow arrangements, (1-3 :2P) and (1-3 :2C), have been graphed 
and fitted by a 5th degree polynomial. Because all computed 
values for effectiveness follow a predictable trend, 
only a sample of the data covering the whole range of 
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values of R have been shown. It is clear that the graphic 
interpretation strongly backs what is known analytically 
from the polynomial regression technique. Where the 
relationship for e = f(N-t u , R) is found explicitly from 

e = A 5 N tu 5 + A 4 N tu 4 + A 3 N tu 3 + A 2 N tu 2 + AjN^ + A 0 (56) 

while the corresponding coefficients A 5 , A 4 , A 3 , A 2 , A^ , and 
Aq relating to a specific value of R are found in Tables 2 
and 3 for the (1-3:2P) and (1-3:2C) configurations. An 
example of how to use this equation in a heat exchanger 
problem now follows. 

B. NUMERICAL EXAMPLE 

Consider a heat exchanger containing 400 m2 of (A = 

400 m^) of surface and operating with an overall heat 
transfer coefficient of 80 W/m2°C. Cold fluid at a capacity 
rate of 10,000 W/°C enters the exchanger at 60°C. Hot fluid 
at a capacity rate of 20,000 W/°C enters the exchanger at 
200 °C . 

1. Find 

The effectiveness ( e) and compute the hot and cold 
fluid outlet temperatures for the (1-3:2C) shell tube pass 
arrangement . 

2 . Assumpt ions 

1) Negligible heat loss to surroundings and kinetic and 
potential energy changes. 
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2) Constant thermal and fluid properties for both 
fluids. 

3 . Analysis 

Here C c = 10,000 W/°C and C h = 20,000 W/°C this 
makes R = C m i n /C max = C c /Ch = (T 4 - T2 )/(t 2 ~ ti) 

10,000/20,000 = 0.5. 

and = UA/C c 

= 80 (400)/ (10, 000) = 3.2 

First, go to Table 2 (page 81) for the (1-3 :2C) arrangement 
with R = 0.5 and find the coefficients 

A 0 = 0.1294 x 10“ 2 

Ai = 0.98120 

A 2 = -0.66161 

A 3 = 0.27938 

A 4 = -0.66456 x 10 _1 

A 5 = 0.66069 x 10 -2 

Then apply equation (56) for N^- u = 3.2 

e = A 5 N tu 5 + A 4 N tu 4 + A 3 N tu 3 + A 2 N tu 2 + A l N tu + A 0 ( 56 ) 
e = 0.769 
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Because C c < 




and with Ti - ti = 200 - 60 = 140°C 

t 2 - ti = e(Ti - ti) 

= 0.769 (140) 

= 107. 7°C 

Finally, the outlet cold fluid temperature t2 is 

t2 = 107.7 + ti 
= 107.7 + 60 
= 167. 7°C 

and the fluid temperature, T 2 is easily found 




T 2 = Ti - 0.5 (t 2 - ti) 

= 200 - 0.5 (t2 ~ t]_) 

= 146. 1°C 

4 . Observations 

The primary observation made here is that by using 
the 5th order polynomial equation (56) with the appropriate 
coefficients found in Table 2 or 3 , an accurate value for 
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effectiveness can be found thus allowing one to solve for 
many more unknown values or prarameters of the heat 
exchanger (i.e., hot and cold outlet temperatures). 

The other observation that is to be made is that 
when comparing the value for effectiveness computed here 
against the value for a 1-2 or 1-4 exchanger (0.745 and 
0.740 respectively) under the same conditions, one finds 
that there is a significant difference in exchanger 
performance as a function of odd or even tube passes and 
that the 1-3 :2C arrangement has a higher effectiveness than 
either the 1-2 or 1-4 arrangement. From inspection of the 
curves for the 1-3 :2P exchanger at Figure N-l or N-6 with 
R = 0.5 and N tu = 3.2 an approximate value of e = .715 is 
obtained. It is clear that this is also less than that of 
1-3 :2C arrangement. Therefore, it is evident that the 
1-3 :2C exchanger out-performs not only the 1-2 and 1-4 
arrangement but also its counterpart the 1-3 :2P exchanger by 
3.1%, 3.8% and 7.0% respectfully. 
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Figure 5.1 1-3 :2C Data Fit by a 5th Order Polynomial 
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Figure 5.2 1-3:2P Data Fit by a 5th Order Polynomial 
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cn 


CM 


in 


ON 


cn 


in 


iH 


I 


cn 


r\ 


iH 




00 


CM 




NO 


m 


ON 


iH 






m 


CO 


ON 


iH 




NO 


I 


NO 


ON 


m 




iH 


cn 


m 


iH 


cn 


iH 


rH 


iH 


iH 


00 


CO 


ON 


r^ 


00 


O 


H 


o 


00 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


00 


00 


00 


00 


00 


00 


< i 
1 


o 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


1 

i 

i 


rH 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 



1 


O 


o 


00 


00 


iH 


ON 


NO 


cn 


o 


m 


ON 




CM 


m 




ON 


ON 


r^ 


m 


n | 


m 


CM 


CM 




iH 


NO 


cn 






00 


r^ 




ON 


r^. 


cn 




ON 


i— i 




O 1 


NO 


NO 


ON 


ON 


00 




NO 


00 


iH 


o 


ON 


m 


00 


ON 


NO 


NO 


o 




o 


rH i 


<r 


ON 


cn 


cn 


o 


NO 




CM 


rH 




00 


cn 


00 


NO 


o 


o 


cn 


NO 


o 


X 1 


m 


o 


CM 


CM 


rv 


cn 


i— i 


NO 






NO 


m 


cn 




NO 




r-v 




ON 


1 

Q 


i— i 


iH 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


< l 

I 


i 






1 

































i 


iH 


CM 


cn 




m 


NO 




00 


ON 


o 


iH 


CM 


cn 


<r 


m 


NO 




00 


ON 


P4 1 
1 


O 


O 


o 


o 


o 


o 


o 


o 


o 


iH 


iH 


iH 


iH 


iH 


iH 


iH 


iH 


iH 


rH 
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TABLE 2 (cont'd) 

- 3 : 2C COEFFICIENTS 



1 

CM | 


00 


CM 


CM 


rH 








CO 


ON 


00 




ON 


CM 






tn 




VO 


00 




O 1 


o 


CM 


CM 


CO 


rH 


CO 


00 


VO 




00 


<r 


00 


to 


rH 


rH 


o 


VO 


o 


ON 


Mt- 


rH 1 


00 


CM 




ON 


CO 


VO 




VO 




00 


VO 


CM 




rH 




rH 




m 




VO 


X 1 




r>* 


00 


00 


O 


o 


rH 


CO 




VO 


rH 


vO 


VO 


00 


ON 


ON 


ON 


o 


o 


rH 


in! 


co 


CO 


CO 


CO 


vj" 












m 














m 


tn 


m 


< ' 
i 


o 


o 


O 


o 


O 


o 


O 


O 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 



O I 


00 




ON 


o 




CO 


rH 


rH 


to 


mT 


00 


rH 




o 




rH 


o 


CO 




00 


rH 1 


m 


rH 


CO 


m 


00 


ON 


00 


O 


m 


o 


o 


00 




CM 


00 


m 


CO 


CM 


00 


m 


X i 


rH 


CO 




00 


o 




CO 


ON 


VO 


m 


VO 


CM 


VO 


00 


00 




CO 


o 


CM 


o 


1 


O 


rH 


CM 


CM 






m 


VO 




ON 


CM 


ON 


ON 


O 


rH 


rH 


CM 


CO 


CO 




^T, 
< 1 
1 
1 


vf 




<r 






<r 




xJ“ 






m 


vJ“ 




m 


m 


m 


m 


m 


m 


m 


O 

l 


o 

1 


o 

1 


o 

1 


o 

1 


o 


O 

l 


o 

1 


o 

i 


O 

l 


o 

1 


o 

1 


o 

1 


o 

i 


o 

1 


o 

1 


o 

1 


o 

1 


o 

i 


o 

i 



i 

i 


O 


ON 


00 


rH 


VO 


o 


O 




CO 


O 


m 


vJ“ 


CM 


CM 


CM 




O 


CM 


CM 




i 


CO 




m 


00 


00 


VO 


rH 


m 


ON 


rH 


rH 


VO 






vO 


VO 


VO 


O 


ON 




i 

ro, 


to 


00 


CM 


CO 




ON 


CM 


VO 


00 




rH 




m 


ON 


CM 


CM 








o 


ON 


ON 


o 


O 


o 


o 


rH 


rH 


rH 


CM 


CO 


CM 


CM 


CM 


CO 


CO 


CO 


CO 


co 




< i 


rH 


rH 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


i 

i 


o 


O 


o 


O 


o 


O 


O 


O 


O 


O 


O 


O 


O 


o 


O 


O 


O 


o 


O 


o 



i 


NO 




rH 


CO 


ON 


rH 


rH 


rH 


o 


rH 




rH 


VO 


o 


00 


CM 




o 


lO 




i 




m 


ON 


CM 


VO 




rH 


Mf- 


m 


rH 


00 




m 


CO 


rH 


VO 


CM 






o 


i 


CO 




CM 




ON 


CM 


m 


O 


CO 


ON 


m 


rH 


CM 




rH 


rH 






00 


CM 


CMi 


to 


m 


VO 


VO 


VO 




r>* 


00 


00 


00 


ON 


ON 


ON 


ON 


O 


o 


o 


o 


o 


rH 


<J l 


m 


m 






m 




m 


m 


m 


m 




m 


m 


m 


VO 


VO 


VO 


VO 


NO 


VO 


1 

l 

l 


o 

i 


o 

i 


o 

1 


o 

1 


o 

i 


o 

i 


o 

i 


o 

1 


o 

i 


o 

i 


o 

i 


o 

i 


o 

i 


o 

i 


O 

i 


o 

1 


o 

1 


o 

1 


o 

1 


o 

1 



1 


CO 


o 


CM 


rH 


CM 


CO 




rH 


CO 


MT 


m 


o 


VO 


<r 


CO 


CO 


VO 


VO 


■4- 


m 


1 


O 


o 


VO 


rH 


ON 


m 


Ml" 


ON 


m 




CM 


00 


CM 


VO 




00 


CO 


CO 


00 


m 


1 








VO 


VO 


VO 


m 


m 


m 




VO 


CO 


CM 


CM 


CM 


o 


o 


o 


00 


00 


rHl 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 






<J ' 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


1 

1 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 



i 

i 


o 


CO 


CM 






ON 


O 




CM 




00 


m 


ON 




O 


CM 


00 


o 


o 


O 


<•» i 




ON 


O 


rH 


CM 


rH 


00 


o 


CM 


CO 




VO 


rH 


CM 


00 


vO 


CO 




CO 




O i 


rH 


00 




m 




rH 


rH 


vO 


r^ 




rH 


00 


CM 


CM 


i— 1 


CM 


O 


ON 


00 




rH 1 


rH 




VO 


rH 


rH 




ON 


ON 


CM 


ON 






CM 


CM 


00 


O 


ON 


rH 


CO 


rH 


X 1 




VO 






VO 


LO 




VO 


vO 




VO 




ON 


00 


00 


rH 


ON 


ON 


o 


rH 


oi 


O 


o 


o 


O 


O 


O 


o 


O 


O 


o 


o 


o 


O 


O 


o 


O 


o 


O 


rH 


rH 



<5 



1 

1 


o 


rH 


CM 


CO 






vO 




00 


ON 


O 


rH 


CM 


CO 


<r 




VO 




00 


ON 


(X 1 
1 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CO 


CO 


CO 


CO 


co 


CO 


CO 


CO 


co 


CO 
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TABLE 2 (cont'd) 

- 3 : 2C COEFFICIENTS 






1 

CM 1 


NO 


1— 1 


CM 


hv 


on 




on 


r>* 


o 


CM 


ON 


00 


on 


NO 


on 


m 


on 




<*■ 


O i 


NO 




NO 


i— 1 


<t 


ON 






on 




NO 


-CI- 


tn 


i — i 


m 


o 


o 




<*■ 


i— 1 1 


m 




CM 


o 


tn 


NO 


on 


on 


T— 1 


i — l 


o 


ON 


tn 


vt 


CO 


<r 




NO 




X i 


i— 1 1 


CM 


on 




on 


on 


tn 


tn 


tn 


O 


NO 


NO 


CO 


00 


ON 


i— i 


o 


i — i 


CM 


i 


to 


m 


tn 


tn 


tn 


tn 


tn 


tn 


m 


ON 


NO 


tn 


m 


m 


tn 


NO 


NO 


NO 


NO 


in 


• 






































< ' 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


O 


O 


o 


o 


o 



i 

H | 

O 1 


O 


00 


00 


o 


00 


ON 


on 




i— i 


NO 


NO 






00 


NO 


1—1 


CM 




o 


1“ 1 1 


on 


CM 


1—1 


ON 


CM 


CM 




o 


CO 


r^> 


m 


ON 


NO 


l—l 


NO 


m 


00 


o 


o 


X i 


i— 1 


00 


NO 


CM 


o 


on 


NO 




NO 


i—i 




on 


NO 




00 


i—i 


<t 


tn 


CM 


i 




<j- 


m 


NO 


NO 


NO 








ON 


NO 


ON 


o 


o 


1—1 


on 


CM 


on 


<r 


Ml 


tn 


tn 


m 


m 


m 


m 


m 


m 


m 


00 


NO 


m 


NO 


NO 


NO 


NO 


NO 


NO 


NO 


< l 

l 

l 


o 

i 


o 

i 


O 

i 


o 


o 

i 


O 

i 


o 


O 

i 


O 

i 


O 

i 


o 


o 

i 


o 

i 


o 

1 


o 


o 

i 


O 

I 


O 

i 


o 

1 



1 

1 


m 


O 


NO 


NO 


00 


NO 


m 


00 


i—i 


r-^ 


CO 


CM 


NO 




CO 




m 


m 


on 


1 


CM 


on 


o 


CM 


ON 


-d“ 


<t 


NO 


CM 


NO 


on 


1—1 




'd - 


r^. 


NO 


CM 


CM 


m 


1 


i—i 


on 


NO 


00 




ON 


on 


on 


-cl- 


00 


ON 


o 


on 




r-^ 




O 


on 


m 




-d- 


-d- 


<t 




-d“ 


m 


m 


in 


m 




NO 


NO 


NO 


NO 


r-. 






r-. 


< i 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


on 


cm 


CM 


CM 


cm 


CM 


CM 


CM 


CM 


CM 


i 


O 


o 


o 


O 


o 


o 


o 


o 


o 


o 


O 


O 


O 


o 


o 


o 


O 


o 


o 



i 

i 


CM 


NO 


on 


CO 


NO 


on 


CM 


m 


NO 


ON 


i—i 


CM 


CO 


CO 


CM 


NO 


i—i 


CM 


CM 


i 


r-. 


o 


ON 


r-^ 


ON 


on 


CM 


<t 




m 


NO 


r-s 


ON 


on 


CM 


ON 


CO 


on 


on 


i 


on 


’-O 


ON 


CM 


CM 


m 


O 


o 


i—i 


i—i 


1—1 


ON 


on 


m 


ON 


on 


CM 


NO 


ON 




i—i 


1—1 


1—1 


CM 


CM 


CM 


on 


on 


on 


r-. 


NO 


on 


-d- 


<1- 




m 


m 


m 


m 


< i 


NO 


NO 


NO 


NO 


NO 


NO 


NO 


NO 


NO 




NO 


NO 


NO 


NO 


NO 


NO 


NO 


NO 


NO 


i 

i 


o 


o 

1 


o 

1 


o 


O 


o 

i 


o 


o 


o 

i 


o 

i 


o 


o 


o 

1 


o 

i 


o 

i 


o 

i 


o 

i 


O 

i 


o 

i 



i 


i-i 


CO 


m 




00 


ON 


NO 


o 


i—i 


1—1 


CM 


NO 


ON 


NO 


1—1 


on 


on 


m 


i—i 


i 




CO 


o 


m 


o 




ON 


i—i 


CM 




O 






CO 




i—i 


o 


r^. 


m 


i 




NO 


r-* 


NO 


m 


-d - 




on 


CM 


NO 


i—1 


o 


o 


ON 


ON 


o 


00 


r-» 


r-* 


i — k 




r-^ 


r-. 


r-» 






r-v 


r-v 




-d - 


00 






NO 


NO 


r-. 


NO 


NO 


NO 


< i 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


i 

i 


o 


o 


o 


O 


O 


o 


O 


O 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 



1 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


O 


o 


o 


o 


o 


o 


n | 


-d- 


CM 


NO 


ON 




i—i 


m 




1—1 


NO 


o 


CM 




t — 1 


NO 


on 


NO 


m 


o 


O 1 




00 


m 


on 


r-^ 




CM 




tn 


1—1 




ON 


i—i 






<1- 


CO 


m 


on 


1 — 1 1 




I— 1 


on 




m 


o 


o 


ON 


m 


on 


ON 


NO 


o 




NO 


NO 


m 


00 


on 


X 1 


r-i 


CM 


i—i 


i—i 


CM 


on 


CM 


on 




NO 


CM 




m 


m 


m 


>d- 


NO 


NO 


NO 


1 

Q 


i—i 


i—1 


i—i 


i—i 


i—i 


i—i 


1—1 


i—i 


i—i 


CM 


i—1 


1—1 


i— 1 


i—i 


i—i 


i—i 


i-i 


i—i 


r— 1 


< ' 
I 









































1 


o 


1—1 


CM 


on 




tn 


NO 




00 


ON 


o 


i—i 


CM 


on 


<3- 


m 


NO 


r^. 


oo 


P* 1 
1 
















<1- 






m 


m 


m 


m 


m 


m 


m 


m 


m 
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TABLE 2 (cont'd) 

- 3 : 2C COEFFICIENTS 



1 

tN | 


o 


vO 


O 


O 


VO 


ON 


4* 


ON 




O 


ON 


00 


o 


CO 


r— 1 


r— 1 


to 


vo 


CO 


ON 


O i 


co 


vD 


o 


r— 1 


ON 


4* 




ON 


CNj 


r— 1 


CO 


r— 1 


4 


CO 


VO 




4 


O 




O 


T— 1 | 


vO 


00 


r^s 


ON 


4- 


VO 


r— 1 


T— 1 


4 


i— 1 


to 


VO 


ON 


to 


4 


ON 


ON 


T— 1 


CO 


4 


X 1 


CN] 


r— 1 


CO 


4 


tO 


VO 






00 


ON 


O 


o 


O 


r— 1 


CNJ 


CN] 


CO 


4 


4 


vO 


1 

LT)i 


vO 


VO 


VO 


VO 


VO 


VO 


VO 


VO 


VO 


VO 


o* 




rv. 














r^s 


< ' 
i 


o 


O 


O 


O 


O 


O 


o 


O 


O 


o 


o 


O 


O 


O 


O 


O 


O 


O 


O 


O 



O i 


to 


o 


4 


vO 


to 


CO 


to 


4 


vO 


4 


4 


CO 


4 


to 


o 


o 


CNJ 


o 


CNJ 


Oh 


T— 1 | 


CO 


CNj 


vO 




o* 


vo 


to 


CO 


CO 




4 


Oh 


CO 


o 


00 


CNJ 


VO 


00 


Oh 


Oh 


X 1 


4 


f— 1 


4 


4 


ON 


ON 


4 


to 


to 


r— 1 


co 


4 


CO 


4 


r— 1 


VO 


4 


vO 


ON 


to 


1 


4 


4 


tO 


VO 


vO 




CO 


CO 


ON 


o 


r—l 


T — 1 


T— 1 


CNJ 


CO 


CO 


4 


4 


"4* 


VO 


^1 
< 1 
1 
1 


vo 


vO 


vO 


vO 


vD 


vo 


vO 


vo 


vO 


Oh 


Oh 


Oh 


Oh 








Oh 


Oh 


Oh 


Oh 


O 

l 


O 

l 


o 

i 


O 

i 


O 

l 


O 

i 


O 

1 


O 

I 


O 

i 


o 

1 


o 

1 


O 

1 


o 

1 


o 

i 


o 

i 


o 

1 


O 

1 


O 

1 


O 

1 


O 

1 



1 


T— 1 


ON 


4 


ON 


00 


VO 


o 


CO 


O 


00 


O 


CNJ 


to 


tO 


r—l 


o 


00 


4 


to 


r>H 


1 


to 


4 


r— 1 


T— 1 


Oh 


00 


vo 


CO 


co 


r— 1 


vo 


CNJ 


VO 


vO 


o 


4 


O 


o 


r—l 


vo 


1 


vO 


vO 


o 


CO 


*4* 




ON 


o 


CO 


to 


00 


ON 


O 


CNJ 


to 


vo 


ON 


o 


i— 1 


to 


C0| 


Oh 


Oh 


00 


00 


00 


CO 


00 


ON 


ON 


ON 


ON 


ON 


O 


O 


o 


O 


O 


r—l 


r—l 


r—l 


< 1 


CnJ 


CNJ 


CnJ 


CnJ 


CNJ 


CnJ 


CN] 


CNJ 


Cn] 


CnJ 


CNJ 


CNJ 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


1 

1 

1 


o 


o 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


o 


o 


O 


o 


O 


o 


O 



1 


00 


4 


r>H 


vO 


tO 


CO 


o 


00 


r—l 


VO 


CO 


tO 


r—l 


vO 


CO 


4 


o 


CNJ 


o 


VO 


1 


vO 


CNJ 


vo 


CO 


CO 


CNJ 


to 


h*H 


co 


VO 


vO 


to 


VO 


CO 


CNJ 


00 


CO 


h*H 


CNJ 


CO 


1 


O 


r—l 


to 


ON 


r—l 


to 


r^H 


00 


CNJ 


4 


00 


ON 


r—l 


4 




CO 


CNJ 


CO 


tO 


o 


CNJ| 


vo 


VO 


VO 


vO 


r^H 




r^H 




00 


00 


CO 


00 


ON 


ON 


ON 


ON 


O 


o 


O 


r—l 


< 1 


vo 


VO 


VO 


VO 


VO 


VO 


vo 


vo 


VO 


VO 


VO 


VO 


VO 


vO 


VO 


VO 




r^H 






1 

1 

1 


O 

i 


O 

1 


O 

I 


O 

i 


o 

i 


O 

I 


O 

i 


o 

1 


O 

i 


O 

1 


O 

1 


O 

I 


O 

i 


O 


o 

1 


O 

l 


o 

i 


o 

i 


o 

i 


o 

i 


1 

1 


VO 


VO 


CO 


to 


to 


r—l 


co 


r—l 


co 


r—l 


r"H 


ON 


r>H 


00 


f^H 


CNl 


vo 


vo 


o 


4 


1 


co 


VO 


4 


co 


to 


r>s 


ON 


r—l 


ON 


CO 


CnJ 


O 


4 


r—l 


CO 


CO 


CO 


ON 


o 


tO 


1 


VO 


4 


to 


to 


4 


4 


CO 


CO 


CNJ 


CNJ 


CNj 


r—l 


O 


o 


ON 


CO 


00 


r*H 


r^H 


p>h 


rHl 


VO 


VO 


vo 


VO 


vo 


VO 


VO 


VO 


VD 


VO 


VO 


VO 


VO 


vo 


tO 


to 


to 


to 


to 


to 


< > 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


1 

1 


O 


o 


O 


O 


O 


o 


o 


O 


O 


O 


O 


O 


o 


O 


o 


o 


O 


O 


o 


o 



1 

1 


o 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


o 


o 


O 


o 


O 


o 


o 


O 


tn | 


r—l 


r^H 


00 


CNJ 


00 


ON 


O 


4 


r—l 


00 




to 


o 


r—l 


ON 


r—l 


4“ 


CNJ 


co 


ON 


O 1 


O 


to 


4“ 


VO 


00 


CNJ 


VO 


CNJ 


to 


00 


CO 


00 


r^H 


CNJ 


VO 


CNJ 


r—l 


CO 


o 


o 


1— 1 1 




to 


vo 


4 


O 


CNl 


tO 


r—l 


CNJ 


r^H 


00 


VO 


00 


O 


to 


00 


O 




4 


o 


X i 




00 


r^H 




00 




00 


00 


00 


00 


00 


ON 


ON 


O 


ON 


O 


o 


o 


r—l 


r—l 


OI 


r—l 


i— i 


r—l 


r—l 


r—l 


r—l 


r—l 


r—l 


r— 1 


r—l 


r—l 


r—l 


1 — 1 


CNJ 


r—l 


CNJ 


CNJ 


CnJ 


CNJ 


CNJ 


< 1 
1 










































1 

1 


ON 


o 


r—l 


CNJ 


CO 


4“ 


tO 


VO 




00 


ON 


O 


r—l 


CNJ 


CO 


4 


tO 


vo 


f^H 


00 


pH 1 
1 


tO 


vo 


VO 


VO 


VO 


vo 


VO 


VO 


vo 


VO 


VO 






f^H 




r^. 




r-H 


r^H 


r^H 



82 



TABLE 2 (cont'd) 

- 3 : 2C COEFFICIENTS 



1 

N | 


CN 


co 


i— 1 




vO 


o 


VO 


On 




CO 


o 


CM 




CO 


CO 


00 


uo 


CO 


vO 


O I 


o 


co 


vO 




i— 1 


vO 




00 


cv 


^d- 


vO 


i— 1 


CO 


CO 


On 


00 


^d- 


O 


On 


i— 1 | 


<r 


00 


CNj 




Mt* 


i—l 


uo 


CM 


CO 


CN 


GN 




CV 


CM 


CM 


vO 


CN 


uo 


CO 


X 1 


VD 




r^. 


00 


00 


o 


o 


r— 1 


f— 1 


i—l 


r—1 


CM 


CO 


MT 


uo 


uo 


uo 


vO 




1 












00 


00 


00 


00 


00 


00 


oo 


00 


00 


00 


00 


00 


CO 


00 


LOl 








































< 1 
1 


O 


o 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


o 


O 


O 


O 


O 


1 

H | 

O l 


i— 1 


i— i 


Cv 


CO 




CO 


CO 


<1- 


vO 


Mf- 


O 


00 


VO 


VO 


uo 


UO 


O 


vO 


uo 


H 1 


GN 


i— 1 


CO 


I— 1 


On 


uo 


i— 1 


CV 


CO 


UO 


^d- 


VO 


r>. 




00 


CO 


CM 


00 


co 


X 1 


vO 


00 


Mf* 


uo 




00 


CO 


00 


o 


uo 


vO 


CM 


CM 


U0 


<r 


00 


i— 1 


U0 


CO 


1 


VO 






00 


00 


CN 


o 


o 


1—1 


r— 1 


f— 1 


CM 


CO 


CO 


<r 


Ml- 


uo 


uo 


VO 


mm 








r>» 


r^. 


r>. 


00 


CO 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


< ■ 
i 
i 


o 

1 


O 

i 


O 

i 


O 

i 


O 

i 


O 

i 


o 

1 


o 

1 


O 

l 


o 

1 


o 

1 


O 

l 


o 

1 


O 

1 


o 

i 


O 

1 


o 

r 


o 

1 


o 

1 



1 

1 


00 


U0 


I— 1 


^d* 


U0 


O 


00 


o 


CO 




I—l 


o 


GN 


r^. 


CO 


O 


CM 


uo 


CM 


1 


CO 


U0 


00 


CM 


CO 


i— 1 


00 


vO 


CM 


ON 


uo 


<1- 


co 


<r 


co 


<1- 


uo 


On 


CM 


1 


o 


ON 


00 


CM 


CM 


vO 




ON 


O 


f— 1 


CM 


<r 




00 


»—l 


CM 


CO 


^d- 


r>. 


COI 


1— 1 


I—l 


f— 1 


CM 


CM 


CM 


CM 


CM 


CO 


CO 


co 


co 


co 


CO 


<r 






<r 


<r 


< 1 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


co 


co 


co 


co 


co 


CO 


CO 


CO 


CO 


CO 


co 


1 

1 

1 


o 


o 


o 


o 


O 


O 


O 


O 


O 


o 


O 


o 


o 


o 


o 


o 


O 


o 


O 



i 


On 


o 


CM 


CM 


CM 


co 


^d- 


^d- 


On 


ON 


vD 


i—i 


VD 


i—i 


vD 


VD 


<r 


VD 


<r 


i 


>d- 


o 




00 


O 


i—l 


vD 


VO 


<r 


vO 


uo 




CM 


vD 


CM 


CO 


ON 


U0 


cM 


i 


i—i 


uo 


<r 


00 


ON 


CO 


uo 




00 


o 


i—l 


CO 




00 


CM 


CO 


<r 


vD 


On 


CM, 


i—i 


i—l 


i—i 


I—l 


i—l 


CM 


CM 


cm 


CM 


CO 


co 


co 


co 


CO 


<r 


<r 


<r 


-d- 


>d- 


< 1 


r>» 








r^. 


r^. 








r^. 


r^. 


r>. 






r>. 


r^. 








1 

1 

1 


O 

i 


O 

i 


o 


O 

i 


O 

i 


o 

1 


o 


O 

i 


O 

i 


o 

1 


0 

1 


O 

i 


O 

i 


O 

i 


O 

i 


o 

i 


o 

r 


O 

i 


o 

i 



i 




CM 


<r 


uo 


CM 


VD 


<r 


CM 


uo 


CO 


CM 


CM 


VD 


co 


VD 


ON 


00 




CO 


i 


<r 


i—l 


UO 


o 


uo 


VD 


CM 


uo 


CO 






O 


O 


o 


CM 


o 


CO 




i—l 


i 


VO 


VD 


Mf 


uo 


co 


CO 


co 


CM 


I—l 


o 


On 


GN 


ON 


00 


00 




vD 


uo 


uo 


i—i i 


uo 


U0 


uo 


uo 


uo 


U0 


uo 


uo 


uo 


uo 


>d- 






>d- 


<r 




<r 


<r 


>d- 


< i 


ON 


On 


On 


On 


ON 


On 


On 


ON 


On 


ON 


ON 


GN 


On 


On 


GN 


On 


On 


ON 


On 


i 

i 

i 


o 


O 


O 


O 


O 


O 


O 


o 


O 


o 


o 


O 


O 


O 


O 


O 


o 


o 


o 



1 


o 


O 


O 


O 


O 


o 


O 


O 


O 


O 


O 


O 


o 


O 


o 


O 


O 


O 


o 






i—l 


uo 


O 


<r 


CO 


On 




i— 1 


CN 


uo 




uo 


00 


vD 


uo 


GV 


<1* 


vD 


O 1 




>d- 


o 




uo 


co 


<r 


ON 


CM 


i— 1 


ON 


CO 


CO 




vD 


o 


<r 


r^. 


uo 


i—l i 


r^. 


CM 


GN 


r^. 


vD 


uo 


i—i 


CO 


CM 


vD 


o 


U0 


00 


CM 


CO 


CM 


co 


o 


CM 


X r 


i—i 


CM 


CM 


i—i 


CO 


CM 


co 


CO 


<1- 


'd" 


uo 


uo 




vD 


uo 


vD 


vD 


r^. 




i 

O 
< i 
1 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


C'l 


cm 



i 


On 


O 


I—l 


CM 


co 


"d" 


uo 


vD 




00 


CV 


o 


r-l 


CM 


CO 


<T 


’U0 


‘sD 




Ul 1 
1 


r^. 


00 


00 


00 


00 


co 


00 


0G 


00 


00 


00 


On 


On 


On 


ON 


On 


ON 


GN 


On 



33 







1 1 
1 M 1 






ON 






i O 1 




m 


o 






1 1— 1 1 


i— i 


o 


i—4 






1 X 1 


CO 


ON 


ON 






i i 


CO 


CO 


CO 






1 LOi 


• 


• 


• 






i <d i 

i i 


o 


o 


O 






i i 

i «-* i 












i O i 


m 




NO 






1 1 — i i 


ON 


CO 


r^ 






1 X 1 


ON 




CO 






1 1 


NO 










1 


CO 


CO 


CO 






l <1 1 


• 


• 


• 






1 l 

1 l 

1 l 


o 

1 


o 

i 


o 

i 






1 1 

1 l 


i—4 


CO 


CM 




c n 


1 1 


CO 




CM 




H 


1 1 


ON 


i—4 


CM 


/— s 


IS 


1 ^ 




m 


m 


T3 


W 


1 <J 1 


CO 


CO 


CO 


► 


H 


1 1 


• 


• 


• 


■P 


CJ 


1 l 


o 


o 


o 


C 


H 


1 1 








O 




l 1 








U 


Pk 


l 1 








V ' 


W 


1 1 


co 


i—4 


NO 




o 


1 1 




CM 


CO 


CM 


CJ 


l 1 




<r 


m 






i <M 


m 


NO 


LO 


W 


CJ 


i <d i 








-1 


CM 


i i 


• 


• 


• 


PQ 


• • 


i i 


o 


o 


o 


<5 


CO 


i i 


i 


i 


i 


H 


1 


i — i 










l—l 


i i 












i i 


CO 


r^ 


NO 






i i 




i—4 


CO 






i i 


<r 


<r 


CO 






i t—i 




<r 








i i 

1 i 


ON 


ON 


ON 






1 1 

1 i 

1 i 


o 


o 


o 



1 


o 


o 


O 


m | 


oo 


oo 


CM 


O 1 


CM 


CO 


CM 


T—I 1 


uo 


ON 




X 1 








1 

o 


CM 


CM 


CM 



<1 



l 

l 


00 


ON 


O 


■ 


ON 


ON 


• 


i 


• 


• 


i—4 
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:2P COEFFICIENTS 



1 

<M 1 VsD 


00 


ON 


i— 1 


NO 


O 


uo 


in 




vO 


m 


r-v 


CM 


Mf 


CM 


i—4 


Mf 


00 


CO 


O 1 NO 


o 


ON 


O 


CO 


m 


O 


i— i 


in 


i— 1 




CM 


Mf 


o 


i— 1 


m 


Mf 


CM 


ON 


i— l i O 




CM 


ON 


CM 


Mf 


CO 


ON 


ON 




i— i 


00 


00 


VO 


CO 


00 


co 


o 


Mf 


X i 


ON 


CO 


Mf 


CO 


00 


CO 






00 


o 


CM 


o 


ON 


1 — 1 


CM 


CO 


m 


m 


1 CM 

lOi i 


i— i 


CM 


cm 


CM 


CM 


CM 


CM 


CM 


CM 


CO 


co 


CO 


CM 


CO 


CO 


CO 


co 


CO 


C i o 


o 


O 


O 


O 


O 


O 


o 


O 


o 


o 


O 


o 


O 


o 


O 


o 


o 


o 



1 

1 1 
O 1 


vO 


in 


co 




Mf 


m 


m 


r— 1 


I— 1 


ON 


o 


CM 


VO 


r— 1 


ON 


CO 


CM 


CM 


i—4 


i—4 I 


in 


Mf 


co 


CM 


VO 


i— i 


NO 


o 


CO 


m 


VO 




VO 


NO 


T— 1 


m 


00 




Ml- 


X 1 


Mf 


CM 


co 


00 


m 


ON 


CM 


m 


f^ 


Mf 


in 


00 


Mf 


m 






00 


CO 


00 


1 


i—4 


m 


00 


ON 


00 


CM 


CO 


CM 


CM 


CO 


Mf 


vO 


in 


Mf 


NO 






ON 


ON 


i 


CO 


CM 


CM 


CM 


CM 


co 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


co 


CO 


CO 


CO 


< i 
i 
i 
i 


O 

l 


O 

i 


o 

1 


o 

1 


O 

i 


O 


O 

l 


o 

1 


O 

l 


o 

1 


o 

1 


o 


o 

1 


o 

1 


O 


o 

i 


o 

1 


o 

1 


o 

1 



1 


NO 


CO 


ON 


NO 


i—4 


o 


CO 


O 


CO 


1—1 


CM 


Mf 


m 


r— 1 


NO 


1—1 


1-^ 


CM 


NO 


1 


Mf 


Mf 


i—4 


CO 


00 


co 


Mf 


CM 


NO 


CO 


C0 


NO 


co 


i-^ 


CM 


m 


ON 




Mf 


1 

CO, 


CO 


Mf 


<r 


ON 


m 


ON 


o 


00 


ON 


CM 


in 


CM 


ON 




CM 


NO 




CM 


Mf 


NO 


Mf 


m 


m 


m 


NO 




NO 


NO 






00 






00 


00 


00 


ON 


ON 


< 1 


rH 


i-^ 


i—4 


r— 1 




I— 1 


i-H 


i—4 


i-^ 


i—4 


i—4 


i—4 


r— 1 


i—4 


1— 1 


1—1 


i— 4 


1—4 


r— 1 


1 

1 

1 


O 


o 


O 


o 


o 


o 


O 


O 


o 


O 


O 


O 


O 


O 


o 


o 


o 


O 


O 



1 




Mf 


ON 


o 


CM 


m 


co 


Mi- 


r—4 


CM 


NO 


CO 


NO 


CM 


ON 


NO 


Mf 


CO 


ON 


1 




o 


CO 


o 


co 


NO 


NO 


oo 


NO 




ON 


NO 


NO 


CM 


CM 


ON 




CM 


n* 


1 


O 


in 


00 


NO 


t— 1 


ON 


o 




O 


Mf 






Mf 


CM 


ON 


Mf 


NO 


co 


m 


CM i 


i—4 


00 


ON 


O 


O 


r—i 


CM 


i-H 


CM 


CM 


CM 


CO 


CO 


CO 


co 


Mi- 


Mi- 


m 


m 


< l 


m 


Mf 


Mf 


m 


m 


m 


m 


m 


m 


m 


m 


m 


m 


m 


m 


in 


in 


in 


m 


l 

1 

1 


o 


o 

1 


o 


o 


o 

i 


o 


o 

i 


o 

i 


o 


o 

i 


o 

i 


o 


o 

i 


o 

i 


o 

i 


o 

i 


o 

i 


o 

i 


o 

i 



i 


O 


m 


ON 


m 






o 


i—4 


CM 


00 


MT 


CO 


1—1 


00 


NO 


Mf 




i—i 


Mf 


i 


CM 


m 


m 


NO 


in 


in 


ON 


CM 


00 


ON 




CM 


Mf 


o 


m 


co 


00 


i—4 


co 


i 


NO 


ON 


Mf 


NO 


r— 1 


ON 


NO 


CO 


CM 


CM 


i-H 


Mf 


l—l 


00 


ON 


o 


00 


O 


ON 


i—l i 


O 


00 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


00 


00 


ON 


00 


ON 


00 


< i 


O 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


i 

i 

i 


i—4 


O 


o 


O 


O 


o 


O 


o 


O 


o 


O 


o 


o 


o 


o 


o 


O 


o 


o 



ft | 


o 


o 


00 


ON 


i—4 


o 




ON 


i—i 


i—i 


1—1 


1—1 


ON 


r— 1 


o 


00 


CO 


00 


CM 


O 1 


m 


00 


m 


CO 


i—4 


m 


ON 


mi- 


i— i 


Mf 








ON 




o 


CO 


o 


CM 


r—4 I 




m 


o 


CO 


CO 


i—4 




nd 




co 


CO 


00 


co 


CM 


00 


o 


00 


Mf 


CO 


X « 


co 


i-^ 


NO 


Mf 


i—l 


NO 


o 


CO 


m 


CM 


CO 


Mf 


O 


o 


m 


ON 


r—4 




Mf 


1 


m 


i-H 


CM 


i— 1 






o 


m 


Mi* 


CO 


NO 


Mf 


Mf 


00 


m 


CO 


NO 


co 


Mf 


oi 
<J i 

i 


r—4 

1 


i — i 


O 


o 


O 


O 

I 


o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


O 


o 


O 



l 


i—i 


CM 


CO 


Mf 


m 


NO 




co 


ON 


o 


i—i 


CM 


CO 


Mf 


m 


NO 




00 


ON 


P4 i 

i 


o 


O 


O 


o 


o 


o 


o 


o 


o 


i—i 


i—i 


i—l 


1—1 


r— 1 


i—i 


l—l 


r— 1 


i—l 


i-H 



85 



TABLE 3 (cont’d) 

-3 : 2P COEFFICIENTS 



1 

N | 


CNj 


o 


CO 


o 


rv 


CM 


rv 


NO 


vr 


vf 


00 


ON 


ON 


rv 


NO 


vf 


r— 1 


Vf 


CO 


rv 


O i 


T—i 


o 


uo 


uo 


ON 


CM 


NO 


NO 


co 


o 


Vf 


o 


CM 


CM 


vf 


NO 


NO 


CO 


00 


CM 


r— 1 1 


o 


CO 


co 


co 


v£> 


Vf 


00 


CO 


CO 


rv 


vX) 


ON 


CM 


o 


NO 


00 


r— 1 


ON 


vf 


co 


X 1 


NO 


rv 


00 


00 


o 


r— 1 


r— 1 


co 


co 


CO 


Vf 


m 


\o 


rv 


NO 


rv 


00 


00 


ON 


ON 


1 

mi 

< i 
i 


CO 


CO 


co 


co 


Vf 


Vf 


Vf 


Ml- 


vf 


Vf 


vf 


vf 


Vf 


vf 


Vf 


vf 


vf 


>d" 


vf 


vf 


O 


o 


o 


o 


o 


o 


O 


o 


o 


O 


o 


o 


O 


O 


o 


o 


o 


o 


o 


o 



O i 


NO 


CO 


o 


rv 


NO 


rv 


uo 


ON 


r—l 


1—1 


vf 


r—l 


CM 


00 


00 


NO 


r—l 


o 


CO 


ON 


1—1 1 


CM 


CM 


ON 


r—l 


rv 


00 


uo 


rv 


r—l 


00 


00 


r—l 


ON 


CM 


NO 


rv 


CM 


co 


uo 


r—l 


X 1 


CO 


vf 


co 


Vf 


CM 


ON 


CO 


m 


IV. 


o 


00 


O 


CM 


O 


00 


00 


CM 


ON 


vr 


Ml" 


1 


O 


r—l 


CM 


CM 


vf 


vf 


uo 


NO 


NO 


rv 


rv 


ON 


ON 


O 


ON 


o 


r—l 


1—1 


CM 


CM 


i 

<C i 

i 

i 

i 


vf 


vf 


vf 


vr 


vf 


Vf 


<1* 




Vf 


vf 


<1- 


Vf 


Ml* 


uo 


vf 


uo 


uo 


NO 


U0 


NO 


o 

i 


O 

i 


o 

i 


o 

i 


o 

1 


o 

1 


o 

i 


o 


o 


o 

1 


o 

i 


o 

1 


o 

1 


o 


o 

1 


o 

i 


o 

1 


o 

i 


o 

1 


O 

i 



I 


00 


uo 


CM 


rv 


uo 


o 


rv 


rv 


ON 


r—l 


r—l 


CM 


o 


CO 


r—l 


00 


co 


NO 


CO 


O 


I 


o 


uo 


00 


o 


Ml- 


o 


i—i 


00 


00 


CM 


00 


uo 


m 


o 


r—l 


r—l 


VO 


o 


ON 


CM 


1 


NO 


ON 


CM 


co 


oo 


r—l 


CM 


m 


NO 


00 


O 


vf 


uo 


00 


00 


r—l 


CM 


uo 


NO 


rv 


fO| 


ON 


ON 


O 


o 


o 


r—l 


r—l 


r—l 


r—l 


r—l 


CM 


CM 


CM 


CM 


CM 


co 


CO 


CO 


co 


CO 


C i 


1—1 


r—l 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


cm 


i 

i 

i 


o 


O 


O 


O 


o 


o 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 



1 


CM 


o 


o 


CM 


vf 


CM 


00 


ON 


00 


00 


vf 


CO 


NO 


uo 


uo 


vf 


00 


ON 


NO 


rv 


i 


00 


vf 


o 


vf 


o 


rv 


r—l 


00 


ON 


NO 


CM 


CO 


NO 


CM 


ON 


rv 


o 


vf 


r—l 


00 


i 


rv 


cm 


rv 


rv 


vf 


rv 


ON 


co 


uo 


rv 


r—l 


NO 


rv 


t-H 


r—l 


uo 


00 


r—l 


Vf 


vf 


CM l 


uo 


NO 


NO 


NO 


rv 


rv 


rv 


00 


00 


00 


ON 


ON 


ON 


o 


o 


o 


o 


r—l 


r—l 


r—l 


< 1 


m 


NO 


uo 


uo 


uo 


uo 


uo 


uo 


uo 


m 


uo 


uo 


m 


NO 


NO 


NO 


NO 


NO 


NO 


NO 


l 

l 

l 


o 

i 


o 

1 


o 


o 

i 


O 

i 


o 

i 


o 

i 


o 


o 

i 


o 

i 


o 

i 


O 

i 


o 

i 


o 


o 

1 


o 


o 

1 


O 

1 


O 


O 

1 


l 

l 


o 


m 


CM 


CO 


ON 


NO 


NO 


rv 


NO 


CM 


NO 


ON 


m 


CO 


rv 


o 


o 


VI- 


CO 


uo 


I 


o 


r—l 


uo 


NO 


NO 


NO 


o 


cM 


NO 


CM 


o 


NO 


CM 


CM 


rv 


vf 


00 


NO 


O 


uo 


I 


00 


00 


00 


NO 


rv 


rv 


NO 


NO 


uo 


vf 


vf 


vi- 


co 


CO 


i—i 


r—i 


o 


o 


o 


00 


i— 1 1 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


rv 


< i 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


I 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 



i 

i 

r> | 


ON 


i—i 


r—l 


ON 


r—l 


uo 


o 


rv 


co 


CO 


CM 


NO 


CM 


CO 


rv 


o 


o 


uo 


o 


o 


O i 


rv 


VO 


UO 


r—l 


00 


uo 


o 


rv 


NO 


CO 


ON 


o 


NO 


ON 


o 


NO 


CM 


ON 


r—l 


NO 


1—1 1 


00 


1—1 


o 


r—l 


CO 


r—l 


CM 


rv 


ON 


r—l 


NO 


CM 


U0 


rv 


uo 


CO 


O 


ON 


00 


rv 


X 1 


uo 


CM 


CM 


O 


uo 


o 


uo 


rv 


VI- 


ON 


CO 


CM 


00 


r—l 


co 


CM 


O 


CO 


CO 


CM 




NO 


NO 


U0 


rv 


U0 


uo 


rv 


NO 


NO 


00 


00 


rv 


00 


00 


ON 


o 


o 


ON 


o 


r—i 


o 


O 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


r—l 


1—1 


o 


r—l 


r—l 



1 


o 


r—l 


CM 


CO 


vf 


uo 


NO 


rv 


00 


ON 


O 


r—l 


CM 


CO 


Vi- 


uo 


NO 


rv 


00 


ON 


pH 1 
1 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


co 


CO 


CO 


CO 


co 


CO 


CO 


co 


CO 


co 
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TABLE 3 (cont'd) 

-3 : 2P COEFFICIENTS 



1 

CM 1 


CM 


CO 


un 


CM 


o 


CM 




VO 


ON 


vo 


o 


un 


ON 


vo 


ON 


on 




CM 


CM 


O i 


ON 




o 


00 


vO 


i — 1 


00 


00 


CM 


rH 






00 






ON 


00 




CM 


rH 1 1 


vO 


cn 


00 


no 


iH 


un 


h 


rH 




on 


on 




ON 






o 


<T 


on 


O 


X 1 


O 


o 


iH 


1— I 


on 


on 




un 


un 


vo 


rH 






ON 


ON 


o 


i—i 


CM 


CM 


1 

LO , 


un 


un 


un 


m 


un 


un 


un 


un 


un 


un 


VO 


un 


un 


un 


un 


vO 


vO 


vo 


vO 


< 1 


o 


o 


o 


o 


O 


o 


o 


o 


O 


O 


O 


o 


o 


o 


o 


o 


O 


O 


O 



i 

H | 

O 1 


vo 


ON 


00 


CM 


00 


ON 


00 




ON 


ON 


ON 


ON 


r^ 


CM 


un 


on 


00 


o 


CM 


1—1 1 


00 


ON 


ON 


i—l 


un 


CM 


on 


un 


NO 


un 




ON 


on 


NO 


r^ 


un 


on 


un 


O 


K i 


VO 


on 


un 


NO 


00 


CM 


ON 




CM 


00 


rH 


00 


<r 


NO 


ON 


on 




CM 


iH 


i 


on 


on 






un 


vo 


NO 




00 


00 


on 


ON 


o 


i—l 


i—i 


CM 


on 






i 


un 


un 


un 


un 


un 


un 


un 


un 


un 


un 


NO 


un 


NO 


NO 


VO 


NO 


NO 


NO 


NO 


< i 


• 






































i 

i 


o 

i 


o 

i 


o 

i 


o 

i 


o 

i 


o 

i 


o 

i 


O 

i 


O 

i 


o 

i 


o 

1 


o 

i 


o 

i 


O 

I 


o 

1 


O 


O 

i 


o 

i 


O 

i 



i 

i 


00 


O 


on 


O 


un 


O 




CM 


un 


NO 




on 


00 


O 


00 


on 


00 


o 


ON 


i 


<r 




on 


00 


NO 


i—l 




CM 


o 


I—l 


on 




CM 


un 


un 


i—i 


1—1 


ON 


ON 


i 


i—i 


o 






00 


o 


cm 


un 




ON 


CM 


CM 


<r 




ON 


rH 




NO 


NO 


00 1 












un 


un 


un 


un 


un 




vo 


vo 


NO 


vo 










<! 1 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


i 


O 


O 


O 


o 


O 


o 


o 


o 


o 


o 


o 


O 


O 


O 


o 


O 


O 


O 


o 



1 

1 


ON 


CM 


r^ 


NO 


ON 






ON 


00 


o 






ON 


o 


i—i 


NO 




CM 


i—l 


1 


NO 


00 


on 


on 


CM 






on 


ON 


ON 


o 




<t 


un 


00 


i-H 


un 


CM 


ON 


1 

CM, 


o 


ON 




un 


O 


CM 


vo 


ON 


iH 




CM 


ON 


CM 




00 


i—l 


<r 


00 


00 


CM 


1—1 


CM 


CM 


on 


on 


on 


on 


<r 


<r 


NO 


<r 


un 


un 


un 


NO 


NO 


NO 


NO 


< i 


NO 


NO 


NO 


NO 


NO 


NO 


NO 


vo 


VO 


ND 


NO 


NO 


NO 


NO 


NO 


NO 


NO 


vo 


NO 


I 

l 

I 


O 

l 


o 

1 


O 

1 


O 

l 


O 

i 


O 

i 


O 

i 


O 

i 


O 

i 


o 

i 


o 

1 


o 

1 


o 

i 


o 

i 


O 

i 


O 

1 


o 

i 


O 

I 


O 

1 



1 


ON 


00 


I—l 


NO 


00 


00 


ON 


rH 


CM 


NO 


ON 


00 


1—1 


rH 


o 


ON 


o 


<r 


iH 


1 


CM 


rH 


<r 


O 




NO 


ON 


cn 


00 


cn 


NO 


ON 


un 


NO 


00 


CM 






un 


1 


ON 






NO 


NO 


un 


un 


un 


<4 


*4- 


o 


CM 


CM 


O 


I—l 


1—1 


o 


o 


ON 


rHl 




r-- 


















00 
















NO 


< 1 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


1 

1 

1 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


o 



c*> | 


O 


o 


O 


o 


o 


O 


O 


o 


O 


O 




o 


o 


O 


O 


O 


O 


O 


o 


O i 


un 




ON 


CM 


o 


CM 


CM 


1— 1 


ON 


00 


cn 


cn 


ON 


un 


CM 


iH 




CM 


cn 


rH 1 


o 


ON 


CM 


O 


1—1 


ON 


O 


I—l 


un 




O 


NO 


cn 


un 


NO 


NO 






CM 


K 1 


i—i 


rH 


cn 


cn 


ON 


un 


vO 


ON 


CM 


vO 


00 


un 


ON 


NO 


un 


ON 


CM 


un 


iH 


1 


i—i 


CM 


i—i 


CM 


rH 


CM 


i— 1 


CM 


CM 


CM 


CO 


cn 


cn 


un 




cn 


un 




un 


o 






































• 


< 1 
1 


i—i 


i—l 


o 


I—l 


i—l 


rH 


i—l 


rH 


rH 


rH 


o 


iH 


iH 


iH 


iH 


iH 


rH 


i— i 


iH 



1 


o 


rH 


CM 


cn 


4 - 


un 


NO 




00 


ON 


o 


iH 


CM 


cn 




un 


NO 




00 


PH 1 
1 










>4 












un 


un 


un 


un 


un 


un 


un 


un 


un 
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TABLE 3 (cont'd) 

- 3 : 2P COEFFICIENTS 



1 

CN | 


o 




00 


00 


H 




o 


On 


ON 


00 


m 


co 


00 


On 


CM 


vO 


CO 


o 


in 


CM 


O 1 


ON 


vO 


i—! 


o 


o 


CM 


rH 


rv 


O 


00 


ON 


co 


CM 


m 


in 


rv 


rH 


vO 


on 


vO 


H 1 


00 


>n 


rH 


ON 


m 




VO 


m 


vO 


I— 1 


ON 


CM 


co 


m 


<t 


i— i 


ON 


ON 


CO 


o 


X 1 


CM 


On 


<f 




in 


rv 


VO 


rv 


00 


ON 


ON 


o 


o 


r— 1 


CM 


CO 


CO 


CO 


m 


VO 


ld! 


vO 


VO 


VO 


VO 


VO 


VO 


VO 


vO 


vO 


vO 


VO 


iv 


IV 


IV 


rv 


rv 


IV 


IV 


iv 


IV 


<3 1 
1 


O 


o 


O 


O 


o 


o 


O 


o 


o 


o 


O 


O 


O 


O 


o 


o 


o 


o 


o 


o 



O I 


VO 


rv 


o 


CO 


ON 


rH 


m 


CO 


m 


o 


o 


CO 


CM 


i— i 


o 


ON 


o 


IV 




rv 


I— 1 1 


CM 


o 


CO 


rH 


in 


00 


rv 




o 


o 




O 


O 


ON 


ON 


rv 


m 


00 


CO 


>0* 


X 1 


00 


O 


o 


iv 


CM 


rv 


CM 


rH 


o 


m 


CM 


m 


IV 


VO 




O 


rv 


00 


o 


VO 


1 


Mf 


O 


vo 


VO 


rv 


00 


00 


ON 


o 


o 


rH 


rH 


rH 


CM 


CO 








VO 


VO 


1 

<3 i 
1 
1 
1 


VO 


rv 


VO 


VO 


VO 


VO 


vo 


VO 


rv 


rv 


rv 


rv 


rv 


rv 


IV 


rv 


rv 


rv 


rv 


rv 


O 

l 


o 


O 

1 


o 


O 

l 


o 

1 


O 

i 


O 


o 


O 

i 


o 

1 


O 

1 


O 

i 


O 

i 


o 


O 

i 


o 

t 


o 

i 


o 

i 


o 

i 



i 


m 


ON 


iH 


m 


CM 


rv 




iH 


co 


O 


rv 


ON 


m 


iH 


CO 


VO 


rH 


ON 


O 


o 


i 


CM 


CM 


m 


VO 


in 


00 


ON 


00 




O 






m 


>0* 


o 


rv 


O 


rv 


iH 


1 — 1 


i 

r0| 


ON 


CO 


CO 


m 


rv 


iH 


O 


CO 


VO 


00 


o 


rH 


cm 


m 


00 


ON 


CM 


CM 


VO 


00 


rv 


ON 


00 


00 


00 


ON 


ON 


ON 


ON 


ON 


o 


O 


o 


o 


o 


o 


rH 


rH 


rH 


iH 


<3 i 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CO 


CO 


CO 


co 


CO 


CO 


CO 


CO 


CO 


CO 


i 

i 


o 


O 


O 


o 


o 


O 


O 


O 


O 


O 


o 


o 


o 


o 


o 


o 


o 


O 


o 


O 



1 


ON 


i — i 


iH 


CO 


CM 


iH 


CO 


ON 


V 


rH 


iH 


IV 


VO 


CM 


VO 


VO 


IV 


'0* 


Mf 


ON 


1 


rv 


o 


VO 


CM 


IV 


00 


rH 


ON 


iH 


iH 


CO 


VO 


vJ- 


IV 


iH 


O 


O 


CO 


rH 


m 


1 


rH 


IV. 


rv 


o 


CM 


IV 


IV 


o 




VO 


ON 


O 


CM 


m 


ON 


iH 




m 


ON 


i— i 


CM | 


rv 


00 


rv 


00 


00 


00 


00 


ON 


ON 


ON 


ON 


O 


O 


o 


O 


rH 


rH 


rH 


rH 


CM 


<3 l 


VO 


VO 


vo 


VO 


VO 


VO 


VO 


VO 


VO 


VO 


vo 


IV 


IV 


IV 


IV 


V 


IV 


IV 


V 


IV 


1 

1 

1 


o 

i 


O 

1 


o 


o 

1 


O 

1 


O 

l 


O 


o 


O 

1 


O 

l 


O 


O 

l 


O 


o 


O 

I 


O 

1 


O 

1 


O 

1 


o 


O 

l 



I 


ON 


MJ* 


o 




m 


m 


m 


"3* 


CO 


IV 


00 


V 


o 


o 


H 


VO 


CM 


CO 


rH 


rv 


I 


o 


O 


CO 


VO 


o 


m 


m 


V 




m 


CM 


iH 


m 


H 


rH 


O 


ON 


ON 


00 


rH 


l 


ON 


CO 


00 


rv 


IV 


IV 


m 


m 


m 




>0* 


CO 


CM 


CM 


CM 


rH 


o 


ON 


ON 


ON 


Hi 


vo 


IV 


VO 


VO 


VO 


vo 


vo 


vo 


vo 


VO 


VO 


VO 


vo 


VO 


VO 


VO 


VO 


m 


m 


m 


< ' 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


i 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


O 


O 


O 


o 


o 


o 



1 

1 

M 1 


o 


O 


o 


O 


O 


o 


O 


O 


O 


O 


o 


O 


o 


O 


O 


O 


O 


O 


o 


o 


o i 


00 


IV 


m 


00 


ON 


00 




00 


O 


rv 


H 


O 


CM 


<r 


IV 


H 


00 


H 


V 


m 


rH l 


in 


vo 


m 


m 


o 


V 


O 


CM 


m 


ON 


in 


on 


O 


VO 


o 


'O’ 


VO 


vO 


V 


00 


X 1 


'O’ 


ON 


00 




ON 


00 


00 


ON 


CM 


VO 


00 


rv 


ON 


CM 


m 


i — I 


rH 


00 


o 


m 


1 

o 
<3 i 


m 


CM 


m 


vo 


VO 


in 


IV 


VO 


rv 


rv 


V 


00 


00 


ON 


00 


O 


ON 


ON 


o 


o 


H 


rH 


H 


H 


rH 


H 


rH 


H 


H 


i — i 


H 


rH 


rH 


H 


H 


CM 


rH 


1 — 1 


cm 


CM 



i 

i 


ON 


o 


rH 


CM 


CO 


<0 


m 


VO 


V 


00 


ON 


o 


i — 1 


CM 


CO 


>0* 


m 


vo 


V 


00 


& i 
i 


m 


VO 


VO 


VO 


VO 


VO 


vo 


VO 


vo 


vO 


VO 


IV 


IV 


IV 


V 




rv 


IV 


V 


V 
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TABLE 3 (cont'd) 

- 3 : 2P COEFFICIENTS 






1 

N 1 




O 




VO 




m 


CNJ 


VO 


CO 


cN 


CO 


<r 


CNJ 


00 


00 


r-l 


CO 


00 


o 


O i 


00 


r— 1 




CN 


VO 


00 




m 


o 




m 


CO 


CNJ 


CNJ 


00 


CN 


00 




CNJ 


r— 4 I 


00 


CN 


00 


m 


m 


m 


i— 1 


r-l 


r-l 


r-l 


i— 4 


00 


O 


r^ 


o 


ON 


CO 


m 


o 


X 1 


LO 


00 




00 


00 


CN 


o 


r—4 


CN 


C\| 


CO 


CN 






VO 


uo 


vO 




00 


1 

LD 1 














00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


<3 i 


O 


O 


o 


o 


o 


o 


o 


O 


O 


o 


O 


o 


O 


o 


o 


o 


o 


o 


o 



1 

H 1 

O i 


CNJ 


CNJ 


r-4 


vO 


00 


r-l 


m 


vO 


O 


ON 


o 




r-l 


<r 


o 


o 


VO 


m 


CO 


i —4 I 


r— 1 


1—4 


00 


ON 


vO 


m 




00 


ON 


o 




ON 


CNJ 


o 


o 


in 


CO 


CNJ 


CO 


X i 


VO 


<r 


CNJ 


00 


ON 


00 


CO 


r-4 


ON 


r-l 


ON 


r^ 


00 




m 


m 


ON 


ON 


CO 


i 


VO 


00 


00 


00 


00 


os 


o 


i—4 


r-l 


CNJ 


CNJ 


CNJ 


CO 




m 


m 


m 


vo 




i 














00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


00 


<3 i 

i 

i 


O 


o 


o 

i 


o 


o 


o 

i 


o 


O 


O 


o 


o 

1 


o 


o 

1 


o 


o 

i 


o 

i 


o 

i 


o 

i 


O 



i 

i 


o 




CNJ 


00 


00 


m 


vQ 


O 


r^ 


CN 


ON 


o 




CNJ 


CNJ 


00 




CO 


<3“ 


i 


<r 


CO 


m 


m 


ON 




CO 


ON 


CO 


O 


vo 


m 


vo 




vo 


r— 1 


m 


m 


00 


i 


00 


CO 


CO 


m 


m 


00 


O 


CNJ 


m 


VO 


00 


00 


r-l 


co 


vo 




00 


r-l 


CNJ 


CO I 


r-4 


CNJ 


CNj 


CN] 


CNJ 


CNJ 


CO 


CO 


CO 


CO 


CO 


CO 


Nf 




nt* 






m 


m 


<3 ' 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


1 

1 

1 


O 


o 


o 


o 


o 


O 


O 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 



i 


r^ 


CO 


CO 


VO 




m 


o 


o 




CN 


00 


m 


CN 


00 


CO 


00 


00 


CN 


r-l 


i 


CN 


vo 


m 


r-l 


VO 


r— 1 


t— 1 


r-4 


o 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Effectiveness values, though not analytically derived, 
can be determined for the 1-3 :2C and 1-3 :2P heat exchangers 
by utilizing a 5th order polynomial approximation. 

Sufficient data now exists to cover a complete range of 
capacity rate ratios for values of N-^u from 0.0 to 3.25. 

From the knowledge of a particular three tube pass heat 
exchanger arrangement and dimensions, fluid flow rates and 
temperature extremes, N^u and R values may be computed. 

Then the effectiveness may be determined using the 
appropriate coefficients from the tables provided herein. 
After the effectiveness is obtained, there are two choices. 

1) From a knowledge of q m ax (see Section III), it is a 
simple matter to determine the actual heat transfer 
rate from the postulation that q = eq max . 

2) Roth fluid outlet temperatures may be computed since 
Q = aj c c ph(t2 “ ti). 

It is also apparent from the work done here that the 
1-3 :2C exchanger effectiveness outperforms that of the 1-2, 
1-4 and its counter part, the 1-3 :2P exchanger, as Ntu 
increases. This is true for any and all values of R. 
Therefore, because it is possible to determine the effec- 
tiveness of a 1-3 exchanger which has a higher effectiveness 
(the 1-3 :2C arrangement) than that of the 1-2, 1-4 and 
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1-3 :2P exchanger particularly at high N^u levels, the 1-3 
exchanger can now be given full consideration in heat 
exchanger design. This would be especially helpful where 
three tube passes could alleviate a configuration problem. 



B . RECOMMENDAT I ONS 

The following recommendations are provided for possible 
follow-on projects of a similar nature. 

• Continue study of the e - N-tu method for the 1-5 heat 
exchangers . 

• Investigate the possibility of using a linear 
approximation from the 1-3 data to find a periodicity 
of R values at which to develop a 1-5 data base. 

• Develop interactive software to be used on a 
microcomputer that contains both the 1-3 and 1-5 
polynomial approximation coefficients. The procedure of 
entering the values for N tu , R, and the type of heat 
exchanger when asked to do so in a menu-driven fashion 
so that effectiveness values can be readily obtained 

is self-evident. 

® After results of 1-5 exchanger analyses are complete, 
investigate the effectiveness for both the 1-3 and 1-5 
heat exchangers experimentally to confirm or refute 
the theoretical work done here. 
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WRITE (IOT, 903) 

903 FORMAT (/ ' INPUT COLD SIDE CAPACITY RATE: 
READ (IN, 902) CCLD 

WRITE (IOT , 904 ) 



904 FORMAT ( / ' INPUT OVERALL HEAT TRANSFER COEFFICIENT 
READ (IN , 902 ) U 
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KCON (1,1) = 1004 
KCON (1,2) = 1014 
KCON (1,3) = 2004 
KCON f 1 , 4 1 = 3015 
COEF ( 1,4) = VALK1 
DO 50 1 = 1,3 
50 COEF (1,1) = VALK3 

NODES 2 TO 50 
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CLOSE ( IWR , IOSTAT=ICK) 



IF(ICK.NE.O) WRITE(IOT, 921)EEX 
921 FORMAT(' Trouble closing printer output file' 
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NTU32P COMPUTER GENERATED INPUT ANALYZER PROGRAM 
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MODIFIED SECTIONS OF THERMAL ANALYZER TO RUN ON BATCH SYSTEM 
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ASSUMING THE USER HAS RUN PROGRAM CHECK AND FIXED ANY WCOEFS 
THAT MIGHT HAVE BEEN IN ERROR. TVSSI, THEREFORE, WILL NOT 
CHECK AGAIN FOR A WCOEF > 100. NO WARNING IS PRINTED. 
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74 WCOEF = WCOEF * ABSTD ** .25 
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MODIFIED NTU14 PROGRAM NTU14BC USED TO RUN ON BATCH SYSTEM 
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APPENDIX M 



1-3 :2C EFFECTIVENESS VS. N fc GRAPHS AT VARIOUS R VALUES 




Figure M.l 1-3 :2C Effectiveness vs. N over Range of 
R from 0.01 to 1.0 



165 




Figure M.2 1-3 :2C Effectiveness vs. N over Range of 

R from 0.01 to 0.10 t 
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Figure M.3 1-3 :2C Effectiveness vs. N over Range of 

R from 0.11 to 0.2 
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Figure M.5 



1-3 :2C Effectiveness vs. N over Range of 
R f rom 0 . 3 1 to 0 . 4 U 
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Figure M.6 



1-3 :2C Effectiveness vs. N over Range of 
R f rora 0.41 to 0.5 
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Figure M. 7 1-3 :2C Effectiveness vs. N over Range of 

R from 0.51 to 0.6 tlJ 
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Figure M.8 



1-3 :2C Effectiveness vs. N ' over Range of 
R from 0.61 to 0.7 
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Figure M.9 



1-3 :2C Effectiveness vs. N over Range of 
R from 0.71 to 0.8 U 
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Figure M.10 



1-3 :2C Effectiveness vs. Ntu over Range of 
R from 0.81 to 0.9 
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Figure M.ll 



1-3:2C Effectiveness vs. N over Range of 
R from 0.9 to 1.0 
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APPENDIX N 



1-3 : 2P EFFECTIVENESS VS. N GRAPHS AT VARIOUS R VALUES 
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Figure N.3 



1-3 :2P Effectiveness vs. 
R from 0.11 to 0.2 



over Range of 
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Figure N.4 1-3 :2P Effectiveness vs. N over Range of 

R from 0.21 to 0.3 
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Figure N.5 1-3 :2P Effectiveness vs. N over Range of 

R from 0.31 to 0.4 U 
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Figure N.6 1-3 :2P Effectiveness vs. N over Range of 

R from 0.41 to 0.5 
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Figure N.7 1-3 :2P Effectiveness vs. N over Range of 

R from 0.51 to 0.6 U 
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Figure N.8 1-3:2P Effectiveness vs. N over Range of 

R f rom 0.61 to 0.7 
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Figure N.9 1-3 :2P Effectiveness vs. N over Range of 

R from 0.71 to 0.8 
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Figure N.10 1-3 :2P Effectiveness vs. N over Range of 

R from 0. 81 to 0.9 tU 
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Figure N.ll 1-3 :2P Effectiveness vs. N over Range of 
R from 0.91 to 1.0 U 
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